Calculation of the anomalous exponents in the rapid-change model of passive scalar 

advection to order e'^. 



L. Ts. Adzhemyan, N. V. Antonov, V. A. Barinov, Yu. S. Kabrits, and A. N. Vasil'ev 
Department of Theoretical Physics, St Petersburg University, Uljanovskaja 1, St Petersburg — Petrodvorez, 198504, Russia 

(24 May 2001) 

The field theoretic renormalization group and operator product expansion are applied to the 
model of a passive scalar advected by the Gaussian velocity field with zero mean and correlation 
function oc 5{t~t') /k'^'^'^ . Inertial-range anomalous exponents, identified with the critical dimensions 
of various scalar and tensor composite operators constructed of the scalar gradients, are calculated 
within the e expansion to order e'^ (three- loop approximation) , including the exponents in anisotropic 
sectors. The main goal of the paper is to give the complete derivation of this third-order result, 
and to present and explain in detail the corresponding calculational techniques. The character 
and convergence properties of the e expansion are discussed; the improved "inverse" e expansion is 
proposed and the comparison with the existing nonperturbative results is given. 
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I. INTRODUCTION 



The investigation of intermittency and anomalous scaling in fully developed turbulence remains essentially an 
open theoretical problem. Both the natural and numerical experiments suggest that the deviation from the classical 
P5 Kolmogorov theory [ 1 1 is even more strongly pronounced for a passively advected scalar field than for the velocity field 
|Jh itself; see, e.g., Ref. [2| and literature cited therein. At the same time, the problem of passive advection appears easier 
tractable theoretically: even simplified models describing the advection by a "synthetic" velocity field with a given 
I Gaussian statistics reproduce many of the anomalous features of genuine turbulent heat or mass transport observed 
' in experiments. Therefore, the problem of passive scalar advection, being of practical importance in itself, may also 
^ [ be viewed as a starting point in studying intermittency and anomalous scaling in the turbulence as a whole. Detailed 
■ review of the recent theoretical research on the passive scalar problem and the bibliography can be found in Ref. [|. 
Most progress has been achieved for Kraichnan's rapid-change model for the first time, the anomalous exponents 
have been derived on the basis of a microscopic model and within controlled approximations [^-^ . 

In that model, the advection of a passive scalar field 9{x) = 0{t,x.) is described by the stochastic equation 



Vt0 = iyod^0 + f, Vt = dt+vA, (1.1) 

where dt = d/dt, di = d/dxi, i^o is the molecular diffusivity coefficient, is the Laplace operator, v{x) = {u,;(x)} is 
the transverse (owing to the incompressibility) velocity field, and / = f{x) is an artificial Gaussian scalar noise with 
zero mean and correlation function 

{f{x)f{x'))=S{t-t')C{r/L), r = x-x'. (1.2) 

The parameter L is an integral scale related to the noise, and C{r/L) is some function finite as L — > oo. 

In the real problem, the field v{x) satisfies the Navier-Stokes equation. In the rapid-change model it obeys a 
Gaussian distribution with zero mean and correlation function 

{v,{x)vj{x')) = Do6{t~t'){2T:)-'^ j dkPy (k) A^^ exp[ik • (x - x')], Afc = e(fc - m) (1.3) 

where Pij (k) = 6ij — kikj /k"^ is the transverse projector, k = |k|, Dq > is an amplitude factor, d is the dimensionality 
of the x space, &{■■■) is the step function and < £ < 2 is a parameter wit h th e real ("Kolmogorov") value e = 4/3. 



The infrared (IR) regularization is provided by the cut-off in the integral (1.3) from below at k — m, where m = 1/i 
is the reciprocal of another integral scale £. The anomalous exponents are independent on the precise form of the IR 
regularization; the sharp cut-off is the most convenient choice from the calculational viewpoints. In what follows, we 
shall not distinguish the two IR scales, setting L ~ The relations 

Do/iyo = .90 = (1.4) 



1 



define the coupling constant go (i-e., the formal expansion parameter in the ordinary perturbation theory) and the 
characteristic ultraviolet (UV) momentum scale A. 

The issue of interest is, in particular, the behavior of the equal-time structure functions 

5„(r) = ([0(t,x)-(?(t,x')r) (1-5) 
in the inertial-convective ran ge A ':$ > 1/r ^ m. 



In the isotropic model (1.1)-(1.3), the odd multipoint correlation functions of the scalar field vanish, while the even 
equal-time functions satisfy linear partial differential equations . The solution for the pair correlation function is 
obtained explicitly; it shows that the structure function 5*2 (x r^~^ is finite at m = [Q. The higher-order correlators 
are not found explicitly, but their inertial-range behavior can be extracted from the analysis of the nontrivial zero 
modes of the corresponding differential operators in the limits £ — > (|J^ ,1/c?— >0|^,ore— >2 . It was shown 
that the even structure functions in the inertial-convective range exhibit anomalous scaling behavior: 

5„(r)oci?o"/'r"(i-^/2)(^^)A„^ r=|x-x'| (1.6) 
with negative anomalous exponents A„ , whose first terms of the expansion in e and 1 /d j^] have the forms 

A„ = ~n{n~2)e/2{d+2) + 0{e^) = -n(n - 2) e/2rf + 0(1/^^). (1.7) 

Another quantity of interest is the local dissipation rate of scalar fluctuations, E{x) = VQdi9{x)diO{x). The equal- 
time correlation functions of its powers in the inertial range have the forms |^J^: 

{E''{x)EP{x')) cx [Kr)-^^^-^^^{mr)^^-^+^^ (1.8) 



with A from Eq. (1.4) and A„ from (pT^). Relations of the form ( |l.^ ) are characteristic of the models with multifractal 
behavior 



In Ref. ||10[ and subsequent papers |10|-|15|1, t he f ield theoretic renormalization group (RG) and operator product 



expansion (OPE) were applied to model (pTl| )-(1.3). In the RG approach, the anomalous scaling for the structure 
functions and various pair correlators is established as a consequence of the existence in the corresponding operator 
product expansions of "dangerous" composite operators, whose negative critical dimensions are identified with the 
anomalous exponents A„. This allows one to construct a systematic perturbation expansion for the anomalous 
exponents, analogous to the well-known e expansion in the R G t heory of critical behavior. 



The key role in the RG and OPE approach to model (L1)-(L3) is played by the critical dimensions A„i, associated 
with the irreducible tensor composite operators 

F^i = inv [d,,e- ■ ■ d,,e {d,ed^0f] , (1.9) 

where I is the number of the free vector indices and n = ^ -I- 2p is the total number of the fields 9 entering the 
operator; the vector indices of the symbol F„; are omitted. The symbol ITZV denotes the irreducible part, obtained 
by subtracting the appropriate expression involving the delta symbols, such that the resulting tensor is traceless with 
respect to any pair of indices. In particular, niV\di9dj9\ — diOdjO — Sij{dkddkO)/d and so on. 



The dimension A„ = A„o of the scalar operator is nothing other than the anomalous exponent in Eq. (1.6); see 



Ref. |l^. The dimensions with I ^ come into play if the forcing (1.2) becomes anisotropic: Ani corresponds to 
the leading zero- mode contribution to the l-th term of the Legendre decomposition for the function Sn', see Ref. [^. 
They can be systematically calculated as series in e: 

oo 

A^.^^A^e'^ (1.10) 



fc=i 



with the first-order coefficient 



n{n-2) {d+l)l{d + l-2) 
2{d+2) ^ 2{d- l)(rf + 2) 



For I = this gives the result of fl, while for n = 3 and Z = 1,3 the results of Refs. |0] are recovered. The result 
(1.11) was rederived later in Refs. [ |16|p7| . 

The coefficients A^^J and were obtained in Ref. |l^ for any n and d; the result for general I is presented in 
H . In particular, one has 
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A|f/ = n{n ~ 2)(0.000203n - 0.02976) - /2(0.01732n + 0.01223) 



(1.12a) 



for d = 2 and 



A^f = n{n ~ 2)(0.00203n - 0.00384) - l{l + l)(0.00710?i - 0.00619) 



(1.12b) 



for d = 3 (analytical results are too cumbersome and will not be given here; see Refs. for / = 0, 2 and for 
general I. 

The 0{e^) contribution to A„; was presented in Rcf. ||l5|| : 

A|f;^ = n{n - 2)(0.00454n2 + 0.06486n + 0.06505) + (--0.01974n2 - 0.10423n + 0.24094 + 0.01748 f ) (1.13a) 
for d — 2 and 

A^f = n{n - 2)(0.00140n2 + 0.01992n + 0.03437) + l{l + 1) [-0.0042071^ - 0.02421n + 0.00280 ^(i + 1) + 0.05065 ] 



(1.13b) 

for d = 3. Here, we have presented the 0{e^) results with improved accuracy in numerical coefficients and corrected 
a misprint in the expression for d = 3 in Ref. ||l^ . No analytical formula for A^^;'' is available for general d, bu t the 
numerical result of the form ( 1.13| ) can be obtained for any given d. The large d limit is discussed in Sec. VIE. 



Besides the calculational efficiency, an important advantage of the RG approach is its universality: it can also be 
applied to the case of finite correlation time or non-Gaussian advecting field; see Refs. [|l^. For passively advected 
vector fields, any calculation of the exponents for higher-order correlations calls for the RG techniques already in the 
0(e) approximation p8|-p0|. Detailed introduction to the RG approach in the statistical theory of turbulence and the 
bibliography can be found in Refs. 22|. 

The main goal of this paper is to give the complete and detailed derivation of the third-order result ( 1.13| ) announced 
in the Rapid Communication |15|, and to present and explain in detail the corresponding calculational technique. It 
might be useful not only for the rapid-change model ( |l.l| )-( pT3| ) and its descendants, but also in a wider context of 
the statistical models of fully developed turbulence and critical dynamics. 

Another scope of the paper is to discuss the nature and convergence properties of the e expansion. The knowledge 
of the three terms allows one to obtain reasonable predictions for finite values of £ ~ 1 and to compare them with the 
existing nonperturbative results: analytical and numerical solutions of the zero-mode equations and numerical 
simulations |2^-^. 

The plan of the paper is as follows. In Sec. |l| we recall the field theoretic formulation of the model, diagrammatic 



technique, ren ormalization and RG equations. In Sec. [II, we briefly discuss the OPE, renormalization of composite 
operators (l.£) and th eir r el atio nship to the issue of anomalous scaling. Since the "ideology" of the RG and OPE 
approach to the model (l.l)-( |l.3| ) is explained in Refs. ]lO|-p^ in detail, here we confine ourselves to only the necessary 



information. In Sec. IV, we present the general scheme of the calculation of the critical dimensions of the operators 
Fni- In Sec. the calculation in the one-loop and two-loop approximations is presented in great detail. Section [V^ is 
devoted to the three-loop calculation; some results for the relevant quantities are given in the Appendix. In Sec. |VIl| 
we discuss the convergence of the e expansion, the improved inverse e expansion and comparison with existing 
nonperturbative results. The main ideas of the paper are briefiy reviewed in the Conclusion. 



II. FIELD THEORETIC FORMULATION OF THE MODEL, DIAGRAMMATIC TECHNIQUE, 

RENORMALIZATION, AND RG EQUATIONS 



The stochastic problem (L1)-(L2) is equivalent to the field theoretic model of the set of three fields $ = {0, 9' , v} 
with action functional 



5($) ^ e'Dge'/2 + 9' [-dt + uod^ - [vd)] 9 



vi?-iv/2. 



(2.1) 



The first four terms in (2T) represent the Martin-Siggia-Rose-type action for the stochastic problem (IT), (1.2) at 
fixed V, and the last term represents the Gaussian averaging over v. Here Dg and are the correlators (1.2) and 
(1.3), respectively, the required integrations over x — (i,x) and summations over the vector indices are understood. 
The model (p.ll) corresponds to a standard Feynman diagrammatic technique with the bare propagators ot, 99, 99' 



(the line 9' 9' is absent). In the diagrams, these propagators are represented by the lines: 
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(2.2) 



the slashed end of a soHd hne corresponds to the field 9' , the end without a slash corre spo nds to 9. The triple vertex 



F($) = ~9'vjdj9 = dj9' ■ Vj9 (the equality holds due to the integration over x in Eq. (2.1)) is represented as 



(2.3) 



The dot with the index j on the solid line denotes the differentiation d/dxj with respect to the argument x of the end 
of the line attached to the vertex; the index j of the derivative is contracted with the index of the end of the line vv 
attache d to the vertex. Owing to the transversality of the vv line, the dot can be moved onto another line, as shown 
in Eq. ( ^.3| ). In the momentum representation the vertex ( |2.3| ) corresponds to the factor —ikj, or, equivalently, +iqj, 
where k and q are the momenta flowing into the vertex via the fields 9' and 9, respectively. The sum of the three 
momenta flowing into the vertex via the flelds 6,9' ,v is equal to zero. 

The line vv in the diagrams corresponds to correlation function ( |l.3| ), and the lines 99' and 99 in model (2^) in the 
a;,k representation correspond to the bare propagators 



{99')^ = {-ico + eu)- 



(99)^ = C(k) {oj' 



(2.4) 



where C(k) is the Fourier transform of the function C from Eq. (1.2) and uj and k "flow via the line from the left to 
the right" if the standard form of the Fourier transform with respect to the coordi nate and time differences is used. 
In what follows, we shall work in the i,k representation, where the propagators ( |2.4[ ) have the forms 

(2.5) 



(99')^ - e{t - t') exp{-(t - t')ek}, {99)^ - {C(k)/2ek} exp{-|t - t'\ e^,]■ 



in {99') t is t he t ime argument of 9 and t' is the argument of 9' . 

The model (2J) is logarithmic (the coupling constant go in Eq. (|l.4| ) is dimensionless) at e = 0, and the UV 
divergences have the form of the poles in e in the correlation functions of the fields 9, 9' . Superficial UV divergences, 
whose removal requires counterterms, are present only in the 1-irrcducible function (9'9^^_.^^, and the corresponding 

counterterm reduces to the form 9'd^9\ see Ref. |l^. Thus for the complete elimination of the UV divergences 
it is sufficient to perform the multiplicative renormalization of the parameters and ~ Do/vq with the only 
independent renormalization constant Z^: 



go = g^i Zg, 



Zr, = Z- 



{Do go^'o = 9^^''t^)■ 



(2.6) 



Here /i is the renormalization mass in the minimal subtraction (MS) scheme, which we always use in what follows, g 
and u are renormalized analogs of the bare parameters go and vo, and Z — Z{g, e, d) a re th e renormalization constants. 
In the MS scheme they have the form "1 + onl y p oles in e." The last relation in ( ^.6[ ) results from the absence of 



renormalization of the contribution with in (2.1 



The renormalized action is obtained from the functional (2.1) by the substitution (2.6) and has the form 



Sr{^) = 9'De9'/2 + 9' [-dt + vZ.d'' ^ {vd)] 9 - vi?-iv/2. 



(2.7) 



where the amplitude Do = g^J-'^v from Eq. (1.3) is expressed in renormalized parameters using Eqs. (2.6). 

The exact response function G = (99'^ satisfies the standard Dyson equation, which in the oj, k representation has 
the form 



G ^(a;,p) = ~iw + j^oP^ - Se/9(a;,p), 



(2.8) 



where the self-energy operator Sg/g in the diagrammatic notation (2.2), (|2.3|) is represented as follows: 



k 

, > - 

'/ \J 

^-H- s \^ 



p-k 



J p 



(2.9) 



The multiloop diagrams, which could be added on the right-hand side of Eq. (2.9), contain effectively closed circuits 



of retarded propagators 



and therefore vanish; it is crucial here that the correlation function {vv) in Eq. (1.3) 
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is proportional to the 5 function in time. Therefore, the self-energy operator is given by the one-loop approximation 
exactly; it is independent of lo and has the form 



Se'e(p) = - 



DoPiPj 



where 



Do (rf-l)p^ 



dkk 



-l~s 



(2.10) 



(2.11) 



with Nk from Eg. (1.3). The quantity 5^ — 2tt''-^^ /T{d/2) is the surface area of the unit sphere in d-dimensional 
space. In Eq. ( 2.1C ) we have used the standard convention Q{t — t') = 1/2 at t — t' for the step function in the 69' 
line and the relation ^Py(k)^ = Sij {d — l)/ d for t he an gular averaging of the transverse projector in d dimensions. 

From the Dyson equation and Eqs. ( 2.10 ) a nd ( 2.11) it follows that the exact response function {96') is obtained 
from its bare counterpart {96')q in Eqs. (2^), (2_^) simply by the replacement 



i/eff = 1^0 + Cdid - 1) Do m ^ /2d£. 



(2.12) 



The renormalization constant Z^, in Eq. (2.6) can be found exactly from the requirement that the "effective diffusivity" 
Veff be UV finite in renormalized theory (2.7), i.e., have no poles in e when expressed in renormalized variables ( ^.(^ ). 
In the MS scheme this gives 



Z^ = \-u{d-l)/2d£, u = gCd 



Cd = Sd/{2n)\ 



(2.13) 



where we have changed to the more convenient coupling constant u. In the renormalized variables, Veg is given by 
the expression 



1 



u{d — 1) {ii/mY 



1 



2d 



(2.14) 



which is obviously finite at e — > 0. 

The basic RG equation for a multiplicatively renormalizable quantity F = ZpFn (correlation function, composite 
operator etc) has the form 



[Drg + if] Fr = 0, Vrg = + - 7.^^?^- 



(2.15) 



Here and below, = xdx for any variable x, and the RG functions (the [3 function and the anomalous dimensions 
7) are defined as 

P = Vf,u, 7F = Vf, In Zp = f3du In Zp for any Zp, (2.16) 

where 2?^ denotes the oper ation at fixed bare parameters go, lyQ. From the definitions, the last relation in ( |2.6| ) 
and exact expression (2.13), for the basic RG functions one obtains: 



/3(u) + 7^], 7^ = u(d- l)/2d. 

From Eq. ( |2.17| ) it follows that the RG equations ( |2.15| ) possess an IR stable positive fixed point: 

= 2de/(d- 1), (3{u^)^0, (3'{u^)=e>0. 



(2.17) 



(2.18) 



This fact implies that correlation functions of model (1.1)-(1.3) in the IR region (Ar ^ 1, nir ^ 1) exhibit scaling 
behavior; the corresponding critical dimensions A[F] = can be calculated as series in e. For the basic fields and 
quantities, including the composite operators 0", the dimensions are found exactly [To| ]: 



(-l + e/2), A[r]=nAfl 



1 -e/2 



(2.19) 



(no c orrect ions o f ord er and higher). This is a consequence of the exact equality 71/ (u*) ~ £, which follows from 
Eqs. (^T7|) and (^T8|) . 

In particular, for the structure functions (1.5) relations ( 2.19| ) along with dimensionality considerations give: 



^„(r) 



(2.20) 



with some nontrivial dependence on the IR scale m contained in the scaling functions ^„(mr). 
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III. COMPOSITE OPERATORS, OPERATOR PRODUCT EXPANSION, AND ANOMALOUS SCALING 



Representations of the form ( 2.20 ) for any scaling functions ^(mr) describe the behavior of the correlation functions 
for Ar ^ 1 and any fixed value of mr. The inertial range corresponds to the additional condition mr <^ 1. The form 
of the functions ^ (mr) at mr ^ is studied using the operator product expansion (OPE) . 

According to the OPE, the behavior of the quantities entering into the right-hand side of Eq. (1.5) for r = x— x' 
and fixed x + x' is given by the infinite sum 



er{t,x) - 9r{t,x.') 



F 



(3.1) 



where Cp are coefficients regular in rri^ and F are all possible renormalized local composite operators allowed by the 
symmetry (more precisely, see below). 



In what follows, it is assumed that the expansion (3.1) is made in irreducible tensors (scalars, vectors and traceless 
tensors); the possible tensor indices of the operators F are contracted with the corresponding indices of the coefficients 
Cp- With no loss of generality, it can also be assumed that the expansion is made in "scaling" operators, i.e., those 
having definite critical dimensions Ap. 

The structure functions ( |l.5D are obtained by averaging Eq. (3.1) with the weight expSa, the mean values (F) 
appear on the right-hand side. Their asymptotic behavior for m — s- is found from the corresponding RG equations 



(2.15) and has the form (F) 



From the RG representation ( 2.20 ) and the operator product expansion (O) we therefore find the following ex- 
pression for the structure functions in the inertial range (Ar 3> 1, mr <C 1): 



5„(r) -i?o""/'r"(i-^/2) 



Ap{mr) {mrY 



(3.2) 



with coefficients Ap regular in m^. 

Some general remarks are now in order. 

O wing to translational invariance, the operators having the form of total derivatives give no contribution to Eq. 
(U): {dF{x)}= d{F{x)) = a X const = 0. 



In model (1.1)-(1.3), the operators with an odd number of fields 9 also have vanishing mean values; their contribu- 
tions vanish along with the odd structure functions themselves (they will be "activated" in the presence of a nonzero 
mixed correlation function (v/ ) or an imposed gradient of the scalar field). 

If the function C in Eq. ( |l.2| ) depends only on r = |r|, the model becomes SO{d) covariant and only the contributions 
of scalar operators enter int o Eq. ( |3.2| ). Indeed, in the isotropic case the tensor indices of the mean values (F) of 
the operators F in Eq. ( |3.l[ ) can only be those of Kronecker delta symbols. It is impossible, however, to construct 
nonzero irreducible (traceless) tensor solely of the delta symbols. 

In the presence of anisotropy, irreducible tensor operators acquire nonzero mean values and their contributions 
appear on the right-hand side of Eq^ ( |3.2[) . In the simplest case of the uniaxial anisotropy, specified by an unit 
vector n in the correlation function ( |l.2| ) , the mean value of a ^ th rank traceless operator is necessarily proportional 
to the Ith. rank symmetrical traceless tensor built of the vector n along with the delta symbols; its contraction 
with the corresponding coefficient Ap gives ris e to the Ith order Legendre polynomial Pi{z) with z — (rn)/r. In 
general, the expansion in irreducible tensors in ( |3.l|) after the averaging leads to the decomposition in the irreducible 
representations of SO{d), the Ith sector corresponds to the contribution of the Ith rank composite operators. 

The leading term in the I th anisotropic sector is given by the I th rank tensor operator with minimal dimension 
A[_F]. The feature typical to the models describing turbulence is the existence of composite operators with negative 
critical dimensions; their contributions in the OPE lead to singular behavior of the scaling functions at mr — > 0, that 
is, to the anomalous scaling. The operators with minimal Ap are those involving the maximal poss ib le n umber of 
fields 9 an d th e minimal possible number of derivatives (at least for small e). Both the prob lem (1.1)-(1.2) and the 
quantities (L5) possess the symmetry 9^9 + const. It then follows that the expansion (3J) involves only operators 
invariant with respect to this shift and therefore built of the gradients of 6. 

In general, the operators entering into the right-hand side of Eq. (3.1) are those which appear in the Taylor 
expansion, and those that admix to them in renormalization. The leading term of the Taylor expansion for Sn is the 
nth rank operator which can symbolically be written as (89)"'; its decomposition in irreducible tensors gives rise to 
operators of lower ranks. In the presence of the noise ( ]l.2| ), operators of the form (89)'' with k < n admix to them in 
renormalization and also appear in the OPE. Owing to the linearity of problem (1.1), operators with k > n (whose 
contributions would be more important) do not admix to the terms of the Taylor expansion for Sn and do not appear 
in the corresponding OPE. 
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We thus conclude that the leading terms of the inertial-range behavior are relate d to the critical dimensions A„i of 
the infinite family of irre duc ible tensor composite operators F„i introduced in Eq. (1.9). 

In general, operators (1.!;) mix in renormalization. One can show that the corresponding infinite renormalization 
matrix 



Fnl = ^ Znl,n'l' 



n'V 



n' V 



(3.3) 



is in fact block-triangular, i.e., Zni^n'i' = for n' > n, and so are the matrix of anomalous dimensions 7(u) = Z~^'D^Z 
and the matrix of critical dimensions A = n + nAg + 7(u*). It is then obvious that the dimensions A„;, given by 
the eigenvalues of the latter matrix, are co mple tely determined by the finite subblocks with n' = n. Therefore, we 

The latter are determined by the 1-irreducible 



al.nl' 



can neglect all the elements of the matrix ( |3.3D other than Z, 
correlation functions with one operator Fni and n fields 9. The diagrams for such functions do not involve the 
propagator (6*^)0 from Eq. (2.4) and can therefore be calculate d di rec tly in the "unforced" model without the noise 
(1.2), that is, without the first terms in the action functionals (2T), (2/7). In the absence of the forcing, the model 
becomes SO{d) covariant, the irreducible operators with different values of I cannot mix in renormalization, so that 
the blocks Zni,ni' appear diagonal. 

We thus conclude that the critical dimensions A„; are determined by the diagonal elements = Zni^ni of the 



matrix (3.3) 



A„i = n + nAe +7„i(u*) = ne/2 + 7„;(u*), 7„i(u) = I?^lnZ„i = /3a„lnZ„; 



(3.4) 



with P from (2.16), from ( 2.1^ ) and Ag from (2.19). Owing to the renormalization, the critical dimension A„; is 
not equal to the simple sum of critical dimensions Ag — (—1 + e/2), Ag = 1 of the fields and derivatives constituting 
Fnl- The elements Zni can be calculated in the model without forcing, in which the operators (1.9) are renormalized 
multiphcatively: F„; = ZniF^i. 



IV. CALCULATION OF THE CRITICAL DIMENSIONS OF OPERATORS F„,: GENERAL SCHEME 



From now on, we shall consider composite operators (1.9) in the model without the noise, that is, with Dg = in 
the action functional (2.1). They are renormalized multiphcatively, Fni = Z^iF^^, and the renormalization constants 
Zni ~ Zni{g, e, d) are determined by the requirement that the 1-irreducible correlation function 



(F«(x)0(xi) . . . 9{xn)),_,^ - Z-l{Fr,iix)9ix,) . . . 9ixn)),_,^ = Z-^T^i{x-x,, . . . ,x„) 



(4.1) 



be UV finite in renormalized theory ( ^.7| ), i.e., have no poles in e when expressed in renormalized variables (2.6). This 
is equivalent to the UV finiteness of the product Z~;^r„;(a;; 0), in which 



Vni{x]9) 



dxi 



dXnTnl{x;Xi, . . .,Xn)9{xi) . ..9{Xn) 



(4.2) 



is a functional of the field 9{x). In the zeroth approximation, the functional (4.2) coincides with the operator Fni{x), 
and in higher orders the kernel r„;(a;; xi, . . . , x„) is given by the sum of diagrams shown in Fig. 1. The analysis of the 
diagrams shows that for any argument Xg, the corresponding spatial derivative is isolated as an external factor from 
each diagram. Using the integration by parts, these derivatives can be moved onto the corresponding fields 9{xs) in 
Eq. (4.2), so that the quantity ([1.2|) can be represented as a functional of the vector field Wi = di9: 



1 



dxi . . 



dx„ rl 



' (x; Xi, . . . , Xn) {Xl) ... Wi^ (x„). 



(4.3) 



The diagrams that determine the kernel F in Eq. (4.3) contain only logarithmic UV divergencies. Therefore, in 



order to find the constant it is sufficient to calculate the functional F with some special choice of its functional 
argument Wi, namely, one can replace it by its value at the fixed point x, the argument of th e op erator F^i in Eqs. ( [1.1| ). 
Now the product Wi-^ {x) . . . Wi^ {x) can be taken outside the integrals over xi, . . . ,Xn in (4.3), so that the functional 
F„i(x; 9) turns to a local composite operator. The integration of the remaining function F„; over xi, . . . , a;„ gives a 
quantity independent of any coordinate variables, and its vector indices can only be those of Kronecker delta symbols. 
Their contraction with the indices of the product Wi-^ (x) . . . (x) gives rise to the original operator Fni {x) with 
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some scalar coefficient F. The integration over Xi, . . . ,Xn means that in the Fourier representation, the corresponding 
correlation function is calculated with all its momenta set equal to zero, which is always implied in what follows. 

Now we turn to the derivation of practical formulas for the calculation of the constants from the diagrams. 
For the sake of brevity, we introduce the notation: 



Zni = Zf, Fni (x) = F, Tni (x; 9)=T = FT. 
Then the UV finiteness of the quantity Z^iTni{x] 9) is expressed by the relation 

Div.p. [Zp^T] =0, 



(4.4) 



(4.5) 



where "Div. p." is the operation which selects the UV divergent part; in the MS scheme, it selects only poles in e. 
Class ifyin g the contributions in model ( ^.7| ) according to the powers of the renormalized coupling constant u from 
Eq. ( p!l3|) gives: 



(4.6) 



k=l 



k=0 



k=l 



where [Z^^]^, F^ and Ffe are the contributions of order u'^ in the respective quantities. We substitute expressions 
(4.6) into Eq. (45) and omit the overall factor F; this gives: 



Div. p. 



Zp^ + Zp^ ^ F,!, 



k=l 



0. 



(4.7) 



We recall that in the MS scheme one has "Z^^^ = 1+ only poles in e," so that Div. p. {Zp^^ = Zp^ ~ 1. Substituting 
this equality into Eq. (17) gives the relation 



Zl 



1 - Div. p. 



Zp^ ^Ffe 



k=\ 



(4.8) 



which allows for the recurrent calculation of the contributions \Zp J in the expansion (4.6) from the quantities F^ 



Indeed, selecting in Eq. (|4.8| ) terms of the same order in u gives: 

{Zp\ = -Div.p. {Fi}, 

[Z^i]^ = -Div.p. {r2+[Z^i]Ji }, 

{Zp\ = -Div.p. {r3 + {Zp\v^ + {ZpXv^\ 



(4.9a) 
(4.9b) 
(4.9c) 



and so on. The relations given in Eq. (19) are sufficient for the three- loop calculation. 



r(i) = 1 




r(2) = 




+ 




No 21 



No 22 



No 31 



No 41 
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r(3) = 




+ 




No 23 



No 32 



No 33 



+ 



No 34 




No 45 No 46 No 51 No 61 

FIG. 1. Diagrammatic representation of the function F in the three- loop approximation. 

In Fig. 1, we present, along with respective symmetry coefficients, all the diagrams needed for the three-loop 
calculation of the function F, except for those with the self-energy insertions of the form (2.9) in the 96' lines. (The 
symmetry coefficients which are not shown are equal to unity.) The index I in F^'-* denotes the number of loops (that 
is, independent integration momenta): ^ = 1, 2, 3. Below each diagram we give its number, which will always be used 
in the following to refer to a diagram. The numbers have the forms "No XY," where X is the number of rays in the 
diagram and Y is the order number (from the left to the right in Fig. 1) of the diagram in the subset with given X. 

The thick dots in the diagrams correspond to the vertices of the composite op erat or F (more precisely, see below), 
all the horizontal dashed lines correspond to correlation functions (vv) from Eq. (F3) with Dq = gh'fi'^, and the "rays" 
correspond to chains of lines {99'){){99')[) ... in the same order on each ray: the upper e nd o f each line corresponds to 
the field 9 and the lower end corresponds to 9' . For this reason, in contrast with Eq. (2.9), in Fig. 1 we do not add 
slashes at the ends of the lines. If desired, they can easily be restored: a slash should be added in the lower end of 
each line belonging to a ray. The lowest "external lines" of the diagrams in Fig. 1 correspond to the factors 9 (and 
not to the propagators). We also note that the diagrams in Fig. 1 do not involve the 99 lines (even in the presence of 
the noise). 

In Fig. 1, w e omitted the diagrams which are topologically possible and wo uld be needed for the general correlation 
function ( |1.3D , but in our model with the delta function in time in Eq. (1.3) they vanish due to the presence of the 
closed circuits of retarded propagators {99') o] cf. the remark in Sec. O below Eq. (2.£). 

The contributions of the diagrams with the self-energy insertions ( |2.9D will autor nati cal ly b e included if the prop- 
agators {99')o in Fig. 1 are t aken to be exact, that is, with Ck = '^effk'^ in Eqs. (2^), (2^). In the renormalized 
variables, Veff is given by Eq. (2.14), the zero order approximation being i/gg = v. It is easy to see that the parameter 
1/ from Ek = I'k^ enters into the final answers for the diagrams as where I is the number of loops in the diagram. 
It is thus sufficient to calculate the diagrams without the self-energy insertions and with ek = i^k^ and then introduce 
the additional factor 



Q 



uid-l)[i^,/mr~^] 

2de 



(4.10) 



for any /-loop diagram; see Eq. (2.14 ). W e stress that the replacement v — > Veff is not needed in the amplitude 
Dq = gvfj,'^ of the correlation function (1.3). 



Expression (2.14) corresponds to the special choice of the function iV^ in Eq. ( |l.3| ). If the specific form of the IR 
rcgularization is different (for example, the function Nk — {k? -\- m^)~'^/^~'^/^_was used in Refs. [5UlO|]), the relation 



( 2.10) for Yjgig remains valid, but the explicit form of the integral U in Eq. ( 2.11 ) changes. Then the quantity Q in 
Eq. (4.10) is given by the following general relation 



Q - [Sfl'fl - Div. p. Sfl'fl] / vp^, 
which recovers the expression ( p^ ) for Nk from Eq. ([J). 



(4.11) 
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I ntro duction of the additional factors (4. IC) to the diagrams of in Fig. 1 with Dq = gvfi'^ in ( |l.3| ) and ek = vk"^ 
in (2^) for the quantities r„ in Eq. (4^) gives 



(1) 



-(2) 



QT 



<1) 



-(3) 



2QT 



<2) 



(4.12) 



Expanding the quantity Q from Eq. (4.10) in e gives rise to contributions with the logarithms ln(/i/m); similar 
contributions also come from the factors (fi/mY^ , which naturally arise in any ^-loop diagram in T^^K It is well known 
that the renormalization constants in the MS scheme are indep ende nt of any m ass p arameters, like m and ^ in the case 
at hand. This means that, after the substitution of relations ( 4.12| ) into Eqs. ( |4.9| ), all contributions with ln(/x/m) in 



Zp will cancel each other. Such cancellation provides a good possibility to control the absence of calculational errors. 

~ [Zp^]^- There, 



In the two-loop calculation jl^ we have checked the cancellation of the contributions s ^ \n{fi/r 



the functi on Nk was taken in the form different from (1.3) (see above), the e expansion of the corresponding quantity 
Q in Eq. ( 4.11 ) contained constant terms along with powers of the logarithms ln{p/m), so that it was necessary to 

take into account the contribution with Q in expression (4.12) for 

Our present choice for in Eq. (1.3) is much more convenient, because the corresponding integral (2.11) contains 



only pole in e, and the e expansion of the corresponding quantity Q in Eq. ( 4.10| ) contains only powers of ln(/i/m) 
with no constant contributions. Since we know in advance that all contributions with ln(/i/m) in Zp^ will cancel each 
other, we ma y simply set /i = m in the calculation of this quantity. Then all the contributions with ln(/i/ m) vanish, 
in Eq. ( 4.10| ) we obtain Q — (no corrections from the self-energy insertions are needed), in Eq. (4.12) we obtain 

= r'' , and the factors (/i/m)''^ in the l-\oop diagrams turn to unity. In what follows, in the calculations of the 

diagrams in F we shall retain the factors (/i/m)'"^ and replace them with unity only in the last step, that is, in the 
calculation of Zp ^ . 

Let us turn to the vertex of the composite operator F, denoted in Fig. 1 by thick dots on the top of the diagrams. 
According to the general rules of the universal diagrammatic technique (see, e.g., Ref. [^), for any composite operator 
F{x) built of the fields 9, the vertex with fc > attached lines corresponds to the vertex factor 



Vkix; xi,...,Xk) = S''Fix)/S9{x,). 



{Xk)- 



(4.13) 



The arguments xi . . . Xk of the quantity ( 4.13 ) are contracted with the arguments of the upper ends of the lines 69' 
attached to the verte x. Fo r our operators (F£), built solely of the gradients Wi{x) = d6{x)/dxi at a single spacetime 



point X, the factors ( 4.13| ) contain the product di^5{x — xi) . . . di^5{x — Xk), and the integrations over xi . . .Xk are 



easily performed: the derivatives move onto the upper ends of the corresponding lines 90' attached to the vertex (such 
derivatives we denote by dots on the lines), and their arguments xi . . . Xk are substituted with x. After the derivatives 
have been moved inside the diagram, the remaining vertex factor for the operator F{x) can be understood as an usual 
derivative: 

V^,...^, (x) = d''Fix)/dw^, (x)... dw,, (x). (4.14) 
In what follows, in order to simplify the notation we shall omit the argument x and use the numerical indices 



, k instead of ii . . . ik- Then the vertex factor (4.14) for a /c-ray diagram takes on the form 

Vu... — did2 ■ ■ ■ F with di = d/dwi. 



(4.15) 



Diagrams Nos 41,46,51,61 in Fig. 1 are "factorizable" in the sense that they can be reduced to the products of 
blocks with a lesser number of rays. All the other diagrams in Fig. 1 will be termed "normal." 



A. Scalarization of the diagrams 



The contribution of a specific diagram into the functional F in (F3) for any composite operator F, built of the 
gradients Wi — diO, is represented in the form 



T = Vu 



Tab. . 
^12.. 



WaWb . 



is the vertex factor (4.15), /"I ' ■ is the "internal block" of the diagram with free indices 



where V12 

WaWb . . . corresponds to external 
number in Eq 



ines. The numerical indices 1,2,. 
(4.16) equals to the number of the letter indices a, 6, . . . 
that is, the number of lines that attach to the vertex of the operator 



(4.16) 
the product 

will always be understood as ji,?2j • ■ -j their 
. and is determined by the number of "rays," 
These lines are given by products of the 
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propagators {O9')o from Eq. (2.5) with ek = vk'^ (see above) and are connected by the hnes vv from Eq. ( |l.3D . As an 
example, we present the general form of a four-ray diagram T: 



1234 




^1234 -^1^34 WaWbWcWd- 



(4.17) 



The thick dot on the top represents the vertex of the composite operator, it corresponds to the vertex factor V1234 in 
Eq. ( 4.15| ). The dots on the lines denote the differential operation di = d/dxf, their indices are shown explicitly on 
t he diagram. The lower external lines of the diagram correspond to factors 9; the derivatives at the external vertices 
(2.3), denoted by dots with the indices a,b,c,d, act on these fields and turn them into the product WaWhWcWd with 
Wi = dO/dxi. After all these differentiations have been performed, all the external momenta in t he diagrams are set 
to zero; the IR regularization is provided by the parameter m in the function {vv) from Eq. (1.3). 

The diagrams are calculated in the time-momentum representation. The integration momenta are assigned to all 
internal lines; the number of independent mom enta is equal to the number of vv lines. The dots with the numerical 
indices 1,2,... on the upper lines in Eq. (1.17) correspond to the vector factors ii/cg, coming from the vertex of the 
composite operator. Here k is the integration momentum flowing via the line, with the index shown near the dot 
(s = 1, 2, 3, 4); the coefficient equals — i if the momentum "flows into the dot" and -l-i if the momentum "flows out." In 



general, similar factors are also present inside the diagram, that is, inside the shaded block in Eq. (4.17). Collecting 
all such factors ±i from the whole diagram gives certain "sign factor" ±1, which, of course, should be taken into 

account in the calculations^ 

Since the vertex factor (4.15) a nd th e product WaWb ■ ■ ■ are symmetrical with respect to any permutations of their 
indices, the quantity " in Eq. ( 4.16| ) is automatically symmetrized with respect to any permutations of the letter 
indices a, 6, . . . and the numerical indices 1,2,.. .. In what follows, such symmetrization will always be denoted by the 
symbol Sym. 

For any fixed number of rays k, the quantity Sym I is represented as a linear combination 



Sym I ^y^ BjSi 



(4.18) 



of certain basis tensor structures Si = {Si)i2"' with certain numerical coefficients Bi. The structures for the fc-ray 
diagrams with k = 2,3,4 have the forms (there are two structures for fc = 2, 3 and three structures for fc = 4) 



k 


= 2 : 


Si 


= Sym [SiaS2b], S2 = Sym [(5i2(Sab]; 




k 


= 3 : 


Si 


= Sym [SiaS2bS3c], S2 = Sym [Si26abS3c]; 




k 


= 4 : 


Si 


= Sym [SiaS2bS3cS4:d], S2 = Sym [Si2SabS3cS4d], 


S3 



S3 = Sym [Si2S346abScd] 



(4.19) 



This is sufficient for the three-loop calculation, because the diagrams with fc = 5, 6 in Fig. 1 factorize into products 
of the blocks with fc = 2,3. 

The quantities which will be directly calculated from the diagrams are not the coefficients Bi themselves, but the 
following scalar quantities related to them: 



A, = tr [{S^rit. Syml^t] = tr [S, ■ Sym I]. 



(4.20) 



Here and below, the symbol tr denotes the contraction with respect to all repeated indices, which will not be shown 

explicitly. 

It is therefore necessa ry to expre ss th e coefficients Bi in Eq. ( 4.1§| ) in terms of the quantities ( 4.2C ). This is easily 
done: substituting Eq. ( fl.lSD into ( [1.20D gives 



Ai ^ M,k Bk , where M,k = tr [S, ■ Sk]. 



(4.21) 



In a compact notation, A = MB (matrix M acts onto vector B and gives vector A). The symmetrical matrix M 
defined in Eq. (4.21) is easily calculated for any given set of structures of the form ( 1.19| ); then the corresponding 
inverse matrix is found, and the desired expressions for B in terms of A follow from the relation B — M^^A. 
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Below we give the explicit expressions for matrix elements Mik = Mki for the structures ( 1.19| ): 



k 


= 2 : 


M 


k 


= 3 : 


M 


k 


= 4 : 


M 






M' 



d{d 
d{d 
d{d- 



l)/2, 

1) (d4 

2) {d^ 



Mi2 = d, M22 = d^; 

2)/6, M12 -rf(d + 2)/3, M22 = d(d + 2)2/9; 
2)(d + 3)/24, M12 = (i(d + 2)(d + 3)/12, M13 = 
l-7d+16)/72, M23 = d(d + 2)2/9, M33 ^ d\d 



d{d + 2)/3, 
-2)2/9. 



(4.22) 



It is worth noting that in the calculation of quantities like Mik in Eq. ( 4.21 ) or Ai in ( 4.2C ), it is sufficient to retain 
the symbol Sym only in one of the cofactors, because the second is symm etrized automatically and can be replaced 
with one of its terms, like the expressions in square brackets in Eq. (4.19). 

Inverting the matrices M in Eq. (4.22) gives the following explicit expressions of the coefficients B in terms of A: 



3 : 



= 4 



Bi = 2a[dAi-A2], 

B2 = a[~2Ai + (d+ 1)^2] with a=[{d- l)d{d - 
Bi =6a[(d + 2)Ai -3^2], 
B2 = 9a[-2Ai + {d+ 1)^2] with a = [{d - 



2r\ 

f 2)(d + 4)]-\ 



l)d{d 

Bi = 24a[{d + 2){d + 4)Ai - 6{d + 2)^2 + 3A3], 

B2 = 72a[-2{d + 2)Ai + {d^ +3d + 6)^2 - {d + 3) A3], 

B3 = 9a[8Ai - 8(rf + 3)^2 + {d + 3)(d + 5)^3] with a = [{d - l)d{d + l){d + 2){d + 4)(d + 6)]"^ 



(4.23a) 
(4.23b) 

(4.23c) 



It should be emphasized that the above relations ( 4.22 ) and ( 4.23 ) are independent of the specific choice of a 
comp osite operator built of the gradients Wi = diO, which only determines the explicit form of the vertex factor 



B. Contractions of basic tensor structures 



Consider the procedure of the contraction of the quantities /fl " in Eq. ( 4.16 ) with external factors: the vertex 
factor V12... of the composite operator and the product WaWi, . . .. General rules given below are valid for any local 
monomial F built of the gradients Wi = diO. The operator F may have vector indices; as a rule, they will not be 
shown explicitly. 

Substituting the decomposition ( [4.18[ ) into Eq. ( 4.16 ) gives: 



T = Y,B, V12... {S^)l^:■ WaWb ■■■=J2 



(4.24) 



where Ti are the contractions of the quantities ( 4.19| ) with the external factors: 

r. = 1^12... {S,)lt. WaWb .... 

Consider first the contractions 

(r,)i2... = {S^Tlt. WaWb ■ ■ ■ 



of the structures Si with the external factors w. For quantities (4.19), they are easily calculated: 



fc = 2 
fc = 3 
fc = 4 



(71)12 = W1W2, (12)12 = w'^Su; 

(ri)i23 = W1W2W3, (r2)i23 = w'^ Sym [612W3]; 

(71)1234 = W1W2W3W4:, (T2)i234 = Sym [S12W3W4], (73)1234 = w"* Sym [612534], 



where the sy mbol S ym d enote s the symmetrization with respect to the numerical indices. 
From Eqs. ( 4.25|) and ( 4.26| ) one obtains: 

r. = i/i2...(T0i2..., 



(4.25) 
(4.26) 



(4.27) 



(4.28) 



where the symbol Sym in the quantities (1.27) can be omitted owing to the fact that the vertices V12... are symmetrical. 
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For a gene ral fc -ray diagram we need to calculate contractions of the vertex factors Vi,,,k with the structures Ti . fc 
of the form ( 4.27 ). Their basis is provided by the set 



Wi 



■ Wk, w'^Si2W3 . . .Wk, w^Si2S34W5 . . .Wk, and SO ou. 



(4.29) 



Consider first the contraction of the vertex factor (4.15) with the first monomial in Eq. (4.29), that is, the expression 
Lk F with the operator 



Lk = Wi . . .Wkdi . . .dk 



di = d/dwi 



(4.30) 



The key role in the calculation of the quantit y Lk F is played by the following consideration: by permutations of the 
factors w and d/dw, the operator Lk in Eq. ( 4.3(]| ) can be represented in the form of a polynomial in the operation 
V = Wid/dwi. The action of the latter onto any monomial F is easily found, namely, "DF — nF, where n is the total 
number of the fields w in the monomial F. By permutations of the factors w and d/dw, one can easily obtain the 
recurrent relation Lk+i — LkiT) — k) for the operator Lk, which along with the relation Li —V gives 

Lk=V(V-l){V~2)...{V-k + l), V = w^d/dwi. (4.31) 

When this operation acts onto F, each symbol V is replaced with the number n, which gives the desired coefficient: 



LkF = n{n - l){n - 2) . . . (n - fc + l)F, 



(4.32) 



where n is the total number of the factors w in the opera tor F. 

Consider now the contraction of the vertex factor ( 4.15| ) with the second structure in Eq. ( 4.29 ). The latter includes 
the factor S12, which after the contraction with 8182 gives the Laplace operator = didi with di = d/dwi. This 
operator commutes with the other derivatives di and acts directly onto F. Therefore the first task is the calculation 
of the quantity d^F, which is easily done for any specific monomial F built of the gradients Wi = did. 

From now on, we shall confine ourselves with the family of operators F^i from Eq. ( |l.9| ), which in the new notation 
take on the form 



Fni = mV [wi 



■ Wl w 



2p] 



n = l + 2p. 



The quantities d^F for such operators are easily calculated: 



d'F^ 



-2,1, 



A„; = {n - l){d + n + l - 2). 



(4.33) 



(4.34) 



It is worth noting that in the calculation of the quantity d^F, there is no need to take into account contributions in 
which the both derivatives act onto the monomial wi . . .wi in Eq. ( 4.33| ): in this case, two factors w are necessarily 
replaced with the symbol Sis (the number of free indices is preserved, and the number of factors w is reduced by 2 
under the action of d^), and any contribution with S jg dis appea rs un der the action of the operation ITZV in Eq. (4.33). 

It is clear that the combination of the relations (4.32) and (4.34) gives a complete solution to the problem of the 
calcu lation of the contractions of vertex factor ( 4.15| ) with basis structures of the form (4.29) for our specific operators 
(4.33), and that the above considerations can be directly generalized to the case of any polynomial in w o pera tor F . 
If the basis structure includes several delta symbols (e.g., two delta symbols in the third monomial in Eq. ( 4.29 )), the 
operation d^ should be applied repeatedly; one should also remember that each next Laplacian acts onto the operator 
which has already been modified (replacement n ^ n — 2 in Eq. ( 4.34 )). After all the operations d^ have been applied, 
only the operation Lk ( 4.30 ) with the reduced number of factors remains; it acts onto the properly modified operator 
F according to the rule ( [4.32 ). For our operators F„;, additional scalar factors uP' , ui*, . . . in monomials ( 4.29 ) restore 
the original scala r fac tor ui^'p in Eq. (4.32). Therefore, all the quantities in Eq. (4.2S) are proportional to the 
original operator (4.32): F^ — kiF, with some numerical coefficients ki. They are easily found using the rules discussed 
above. In particular, for the 2-ray diagrams Fi — V12 W1W2 — L2F — n{n — l)F, F2 — V12 w^Si2 ~ w^d'^F — XniF, 
and similarly for the 3-ray and 4-ray diagrams. 

Substituting the relation F^ = kiF into Eq. (4.24) gives: 



F = ^F, T = Y,hB„ 



(4.35) 



where ki are the coefficients in front of F in the contractions of the vertex factor Vi,,,k with structures (4.29), numbered 
in the same order (fci corresponds to the structure without delta symbols, k2 corresponds to the structure with one 
delta symbol, and so on): 
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(4.36a) 
(4.36b) 
(4.36c) 
\ni; (4.36d) 

(4.36e) 

The coefficients ki with k — 5 and 6 will be needed only for the calculation of the factorizable diagrams Nos 51,61 
with five and six rays, repectively. 

Th e quantities Ai are calculated directly from the diagrams, then one calculates the quantities Bi u sing the rel ation s 
(4.23), and then, finally, one calculates the desired contribution of the diagram T using the relations (4.35) and (4.36). 



k = 


2 


fci 


= n(n 




k2 = Kiu where A„i = [n- L){d + 


n + I — 


2); 


k — 


3 


fci 


— n(n 




-2), fca (n - 2)A„;; 






k = 


4 


fci 


= n(n 


-l)(n 


-2)(n-3), fca = (r^-2)(n-3)A„^ 




An-2,iAni; 


k = 


5 


fci 


= n(n 


-l)(n 


-2)(n-3)(n-4), fca = (?i - 2)(n - 


-3)(n- 


4)A„;, /C3 = (n-4)A„ 


k = 


6 


fci 


= n(n 


-l)(n 


-2)(n-3)(n-4)(n-5), fc2 = (n - 


-2)(n- 


3)(n-4)(n-5)A„,, 






fc3 


= (n — 


4)(n- 


5)A„_2,Aji(, ^4 = A„_4_;A„_2,;Ani- 







C. Calculation of the quantities Ai 



The quantities " for the diagrams T are calculated at zero external momenta; the independent integration 
momenta will always be denoted by k in the one-loop diagrams, k, q in all two-loop diagrams, and k, q, 1 in all 
three-loop diagrams. They are always assigned to the horizontal lines vv; in the normal (not factorizable) diagrams 
the order is the following: k flows via the uppermost line, q flows via the next line, and 1 flows via the lowest line. 

These rules are illustrated by Fig. 2, where the one-loop and two-loop diagrams are presented with all their momenta; 
the directions of the latter arc shown by arrows. We denote all the propagators by solid lines, because the horizontal 
vv lines cannot be confused with the "rays" built of the 69' lines. The letter indices a, 6, ... of the vv lines are always 
free (they will be contracted with the indices of the external factors WaWh . . .), while the numerical indices of the vv 
lines are always contracted with the indices of derivatives (that is, the momentum factors ±ips), shown by dots on 
the 69' fines. 




No 21 No 22 No 31 No 41 



FIG. 2. Diagrams of the function F in the one-loop and two-loop approximations in the detailed notation. 



We sfiall calculate the diagrams in the time-momentum (i-k) representation. First, the elementary integrations 
over all times are performed. This gives certain "energy denominator" (f>E, that is, the factor (f>~^^ in the integrand. 
The remaining integrations are those over d-dimensional momenta (k in one-loop diagrams, k, q in two-loop diagrams, 
and k, q, 1 in three-loop diagrams). Each integration is accompanied by the factor Do{2tt)^'^ with Dq — gi'iJ,'^, and 
all z^'s will later be cancelled out with analogous factors in the energy denominators. As a result, the quantity 1^2 
from Eq. ( 4.16| ) is represented as follows: 



1 loop : 

2 loops : 

3 loops : 



jab. 
^12.. 



Tab. . 
-'12.. 



rab. . 
-'12.. 



= [i?o(2vr)-'*] J dk 
= [Z?o(2^)-'^]' J Jdkdq 
= [L>o(27r)-'^]^ / // '^^'^^^^ 



ah. . . 
12... I 



(4.37a) 
(4.37b) 
(4.37c) 



with from Eq. (1.3). After these scalar factors have been isolated, only the transverse projectors P remain in the 
vv lines. Tfie nontrivial parts (j)'^" of tfic integrands in Eq. (4.37) are represented as follows: 



^12... 



rab. . 
-'12.. 



(4.38) 
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where C is a numerical factor, 0£ is the energy denominator, c^g is an additional "scalar factor," which can arise due 
to contractions over the internal indices, and I12 " is the nontrivial "index structure" of the integrand, which remains 
after (hg has been isolated. 



Let us explain in more detail the calculatio n of the factors entering into Eq. ( 4.3j ), starting with the energy factors 
4>~^^ ■ Since the velocity correlation function (1.3) contains the delta function in time, the dots on the rays connected 
by such correlators have coincident times t. If we assign the time t = to the uppermost vertex of a diagram (that is, 
the vertex of the composite operator), the independent variables will be the times < of the horizontal vv lines. On 
any of the rays, the times ti decrease from above to below owing to the retardation property of the functions {09')q 
in Eq. (2.5). For the normal (not factorizable) one and two- loop diagrams in Fig. 1, the order of the times < is 
determined by the form of the diagram in a unique way, that is, there is only one "temporal version." This is also 
true for all the three-loop diagrams with two and three rays. However, for two of the three-loop diagrams (Nos 43,44 
in Fig. 1) there are two equivalent temporal versions (the exchange ^2 ^3 in diagram No 43 and ti ^ <2 in diagram 
No 44). This gives the additional factor 2 in the coefficient C in Eq. (4.38) for these two diagrams. 

The energy denominator cjiE, obtained as the result of the integrations over all time variables, is given by the 
product of the factors corresponding to the "temporal cross-sections" of the diagrams between the horizontal vv lines; 



to each cross-section corresponds the sum of "energies" 



^ for all intersected 96' lines. For example, one has 



(pE = 2ek for diagram No 21 in Fig. 2 
No 31. 

The scalar factors 



iiE — 2eq • 2ek+q for diagram No 22, and (j)E — 2eq [ek + Cq + fk+q] for diagram 



in Eq. (4.38) are generated by the vv lines whose vector indices are contracted with the indices 
of two integration momenta in the diagram. This gives rise to expressions of the form (pPp') = piPijPp where P is 
the projector that corresponds to the vv line and p, p' are the momenta contracted to the line (the factors ±i are 
not included). In what follows, we denote as Py^, Pq and Pi the transverse projectors for the momenta k, q and 1, 
respectively. Then the scalar factors from the diagrams in Fig. 2 take on the form (j>s = (qPkq) for No 22 and 0s = 1 
for Nos 21,31, because the latter two di agra ms do not contain lines vv contracted to two momenta. 

The numerical coefficient C in Eq. ( 1.3^ ) is represented as the product C = C1C 2C3, where C'l is the symmetry 
coefficient of the diagram, shown in Fig. 1, C2 = ±1 is the sign factor (see Sec. tV A ) and C3 = 1 or 2 is the number 
of the temporal versions. For the most of the diagrams one has C2 = +1, but there are two diagrams with C2 = — 1 
(Nos 44,45). The factor C3 is equal to 1 except for the diagrams Nos 43,44 with C3 = 2. In what follows, we shall 
always present the coefficient C as the product C — C1C2C3, for example, C = 1/2 • 1 • 1 for diagrams Nos 21,22 in 
Fig. 2 and C = 1 • 1 • 1 for No 31. _ 

have 



Now let us turn to the quantities Jfl' " in Eq. (4.38), which remain in the integrands after the scalar factors 



been isolated. They have the same stucture for all the diagrams with two and three rays and any number of loops, 
and two possible structures for the four-ray diagrams, namely: 



fc = 2 

fc = 3 

fc = 4 

fc = 4 



nl = I^^ = X^X2Pab-: 

n2 = = xi{x + y)2V?,aaPbc] 

n3 = /123I = xiy2Zs.{x + v + z)iaal3bPcd] 

n4 = = xiy2Z:i{x + y + z)iPabPcd- 



(4.39) 



We have introduced the numbering nl-n4 and in what follows, when describing the diagrams we shall only give the 
number of the corresponding structure, e.g., I = nl with the specification of the vectors x, y,z,Q;,/3 and projectors 
P, P'. In particular, we have I — nl with x = k, P = Pk for diagram No 21 in Fig. 2, I — nl with x = k + q, P = P^^ 
for diagram No 22, and I = n2 with x = k, y = q, ct = Pkq, P = Pq for diagram No 31 (quantities like a = Pkq are 
undestood as vectors obtained by the action of the matrix Pk o nto the vector q) . 

Our aim is the calculation of th e quantities Ai in Eq. (4.2C) for ea ch di agram in Fig. 1. They are given by the 
contractions of the i ntegr als (4.37) with the structures Si from Eq. (4.19), which leads to the replacement of the 
in Eq. ( 4.37 ) with the quantities = tr [iSi)i2.'.'. ■ Sym 4'i2.'.'.\ ■ Substituting expression ( 1.3^ ) into this 



factors 012!! 



construction gives 0^ = Cff'E where 



Ai = tr [S.i ■ SymT\ 



(4.40) 



are t he an alogs of quantities ( 4.20 ) with the replacement I ^ I. As a result, from the relations ( 4.20| ), ( 4.37 ), ( 4.38| ) 
and (4.40) one obtains the following explicit formulas for the calculation of the coefficients Ai: 



1 loop : A, = [Do{2t:)-''] J dkNkCcj)^^ cj), A, 



2 loops : A, = [Do{2t:)-'^]' 



dkdciNkNqCcj)-^ cj). A, 



(4.41a) 
(4.41b) 
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3 loops 



A, = [Do{2T:)-''f 



dkdqdl NkNgNi C I 



: Ai 



(4.41c) 



with Ai from Eq. (|4.40| ). 

Calculation of these quantities for structures nl and n2 from Eq. ( 4.39| ) and Si from (4.19) gives: 

fc = 2(/==nl): Ai = (xPx), A2 ^ {d - l)x'^; (4.42) 
fc = 3 (/ - n2) : 3Ai = (ya)(xPx) + (xa)(yPy) + ((x + y)a)(xPy), 

9^2= [a;2 + 2(xy)][(d-l)(ya) + 2(yPa)] + [y^ + 2(xy)] [(d ~ l)(xa) + 2(xPa)] . (4.43) 

Analogous expressions can also be obtained for the stru cture s n3 and n4 in Eq. ( 4.39| ) for the four-ray diagrams; 



they are rather cumbersome are will be given later in Sec. VI A 



V. CALCULATION IN THE ONE-LOOP AND TWO-LOOP APPROXIMATIONS 



In the one-loop approximation, we need the only diagram No 21 , and in the two-loop approximation, all the four 
diagrams from Fig. 2 a re nee ded. Using the general rules discussed in the previous Section, for the normal diagrams 
in the notation ( 4.38| ), ( 4.39 ) we obtain: 



No 21 : 


C = 


1/2 


1-1, 


4 


>s = 1, = 2ek, 




/ = 


nl 


with 


X 


= k, P = Pk [Px = 0]; 


No 22 : 


C = 


1/2 


1-1, 


4 


)s = (qPkq), 4>E ^ 4eqek-Hq, 




I = 


nl 


with 


X 


= k + q, P = Pq; 


No 31 : 


c = 


1 • 1 


•1, ' 


Ps 


= 1, 0s = 2eq(ek + Cq + ek+q), 




I = 


n2 


with 


X 


= k, y = q, a = Pkq, P = 



(5.1) 
(5.2) 

[(xa)==0, Py-0]; (5.3) 

in the square brackets we give simple relations specific of a given set of momenta and projectors, which are useful in 
the calculations. We also recall that ep =j/p^ 



Substituting the explicit expressions (|5.l| j-(pT^) into the general relations ( 4.42 ), ( 4.43| ) gives 



No 21 
No 22 
No 31 



ii=0, A2 = (d-l)fc2; 

il = (kPqk), i2 = (d-l)(k + q)2; 

3ii = (qPkq) (kPqk), 9^2 = id-l) (qPkq) [k'' + 2(kq)] + 2(kPqPkq) [q^ + 2(kq)] . 



(5.4) 
(5.5) 
(5.6) 



A. One-loop approximation 



The one-loop approximati on F ^ in Eq. (^^) is determined by the only diagra m No 21 in Fig. 2. Substituting the 
quantities known from Eqs. (5.1), ( |5.4| ) and the relation Ck = vk^ into Eq. (4.41a) gives (we recall that the symmetry 
coefficient 1/2 is included into C): 



No 21: Ai^O, A2 = [(d-l)/A]gfi''U 



(5.7) 
(5.8) 



with the integral U = (2tt) J dkNk from Eq. ( 2.11 ). Substituting that expression into (5/?) gives: 

No 21 : ^1=0, A2=u[{d-l)/4:e]{^i/my, u^gCd 

with coefficient Cd fr om Eq . ( 2.11 ) 

From the relations (4.23a ) and (5.8) we find Bi — ~2aA2, B2 = a{d+ 1)^2 with a from (4.23a) and A2 from ( ^.^ 
then from Eqs. ( 4.35 ) and ( 4.36a ) we obtain the quantity F for our diagram: F(No 21) 0:^2 [— 2n(n— 1) + ((i+ l)A„i 
Substituting the known expressions for a and A2 gives: 

(d+l)A„,-2n(n-l) 



F^^^ = F(No 21) = u{n/my 



Aed{d + 2) 



Then the relation (4.9a) gives the final result for the 0{u) contribution to the renormalization constant: 

2n(n- 1) - {d+l)Xni 



Aed{d + 2) 



(5.9) 



(5.10) 
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B. Two- loop approximation 



Co nsider now the two-loop diagrams Nos 22,31 in Fig. 2; the corresponding coefficients Ai are known from Eqs. 
(U) and (l^l). They include the quantities (kPqk) = Psin'^'d, (qPkq) = (J^ sin^ i9, (kPqPtq) = -(kq)sin2i9, 
(kq)^ ~ fc ^g"^ cos^ 1?, where is the angle between the vectors k and q. Substitutuing these expressions into Eqs. 
( ^.5D and (5^) gives: 



No 22 : 
No 31 



Ai=k'^sm^d, ^2 = (d- l)(k + q) 
3li = fc\^ sin-* I?, 
9A2 = sin^^{{d-l)q^[k 



g2 sin^T^Kd- l)[fc2 



2(kq)]-2(kq)[q2 + 2(kq)]} 
2(kq)] - 2(kq) - 4k^ cos^ i?} 



Substituting these expressions along with the other needed information from Eqs. 
Ep = i^P^ into Eq. (4.41b) gives 



No 22 



A, = [(gA^-)V8](2^)-2'^ 
A2 = [{d-^){g^^'f/8]i27r) 



dkdqNkN, 



fc^ sin** 



-2d 



'(k + q)2' 
dkdciNkNq sin^??; 



(5.11) 



(5.12) 



(5.3) and the relation 



(5.13) 



No 31 



3Ai = [(5M^)V4](2^)-2'^ 
= [(.gM^)V4](2^)-2'^ 



fc^ sin'* 'd 
^^^q^^-^^ P+g^ + (kq) ' 



dkdqNkNr 



sin^i9{(d- l)[fc2 + 2(kq)] - 2(kq) - 4fc2 cos^ i?} 
A:2 + q2 + (kq) 



(5.14) 



We note that in the coefficients Ai,2 for diagram No 31 the factor from eq in the denominators cancels out with 
the analogous factor in the numerators, while in the coefficient A2 for diagram No 22 a similar cancellation of the 
factor (k + q)^ takes place. 

The denominators in the integrals (5.13), ( |3.14 ) are symmetrical with respect to the permutation k ^ q. This 
allows one to perform the analogous symmetrization in the numerators, which is equivalent to the replacement 
P (fc2 + g2)/2 in the latter. Then the expressions simplify and reduce to linear combinations of the following 
standard integrals: 



2p 



(27r) 



-2d 



J J dkd(iNkNq sin^P 1}, P=l,2, 



(5.15) 



(2^) 



-2d 



J J dkdq NkNq (kq) sin'' /(k + q)^ 



(5.16) 



h2p = (27r) 



-2d 



J J dkdqNkNg (kq) sin^" ^ /[P + + (kq)], 



p=l,2. 



The coefHcients A, in Eqs. (|t|), (| TJ) are expressed in terms of these integrals as follows: 
No 22 : A, = [{g^i' f/8] [H,/2 - h], A^ = [{gfi')V8] [{d-l)H2]; 

No 31 : 3Ai = [(g/i")V4] [H4/2 - hi/2], 9A2 = [{gfi'f/i] [{d - b)H2/2 + 2Hi + 3(rf - l)/i2/2 - 2hi] 



(5.17) 



(5.18) 
(5.19) 



These expressions along with Eq. ( 4.23a[) for diagram No 22 and Eq. (4.23b) for diagram No 31 give the coefficients 
Bi, then using Eqs. ( 4.35 ), ( 4.36a ) and ( [4.36b ) one obtains the correspon ding qu antit ies r(No 22) and r(No 31); the 
answers are represented as linear combinations of the standard integrals ( 5.15| )-(5.17). 



-(2) 



The two-loop quantity T also includes the contribution from the factorizable diagram No 41 with the symmetry 
coefficient 1/8; see Fig. 1. It reduces to the contraction of the vertex factor "V1234 with the product T12T34 of two 



identical one-loop blocks similar to diagram No 21. In the notation of Eqs. (4.24)-(4.25 ) this contribution has the 
form 
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(5.20) 



where B12 are the known coefficients for the one-loop diagram No 21 (see Sec. V A) and C is an additional symmetry 
coefficient; in the case at hand, C = 1/2. Let us explain its origin. By d efinit ion, the quantities B12 for diagram 
No 21 already include its symmetry coefficient 1/2 (in the factor C in Eq. ( 4.38 )), which gives 1/4 for the product of 
two such blocks. Therefore, in order to obtain the correct symmetry coefficient 1/8 for diagram No 41, the additional 

coefficient C = 1/2 should be included. 

Multiplying the ex press ions in the square brackets in Eq. (5.20) and taking into account the definition of the 
coefficients ki in Eq. ( 4.36 ) gives 



r(No 41) = [1/2] [kiBf + 2k2BiB2 + k^Bj 



(5.21) 



with coefficients fci.2,3 from Eq. ( 4.36c ) and known quantities Bi,2 for the one-loop diagram No 21. In what follows, the 
same scheme will be used for the cal culati on of t he contributions of the three-loop factorizable diagrams Nos 46,51,61; 
the coefficients ki from Eqs. ( 4.36cj ) and ( 1.36c ) will be involved. 
We have found all terms in the two-loop expression 



r^^^ = r(No 22) 



r(No 31) +r(No 41). 



(5.22) 



Now Eq. (4.9b) is used to determine the two-loop contribution [Zp^]^ in the renormalization constant Zp^, with 



r^^^ from Eq. (5.9), F^^"* from Eq. (5.22) and [-^i?^]]^ from Eq. (5.1C). It is also necessary to set /i = m in all the 



diagrams of ' ' in order to eliminate contributions with the logarithms ln(/i/m) (they would be cancelled with the 
contributions of the diagra ms wit h self-energy insertions, which we already omitted). 

All the quantities in Eq. (4.9b) are known explicitly, except for the contributions of diagrams Nos 22,31, for which 
we only know the expressions of the coefficients Ai in terms of the standard integrals (5.15)-(5.17). In order to 
calculate the quantity (4.9b), we need only their divergent parts ~ 1/e^ and 1/e, but for the three-loop calculation 
we shall also need the zero-order in e (constant) terms. 

Calculation of these standard integrals is a separate task and will be discussed in the next Subsection. 



C. Calculation of the two- loop integrals ( 5.15 )— ( 5.17 ) 



In this Subsection we denote — k/fc for any vector k, diik is the area element of the unit d-dimensional sphere 
and (. . .) is the averaging over the sphere. In particular. 



dkNk 



dk 



diik 



Sd 



k 



dk , ^ 
T+7 V ■ ■/ 



(5.23) 



with Nk from Eq. (1.2); Sd — 2n'^^^ /T{d/2) is the surface area of the unit sphere. 

For any two vectors with the angle between them, the following formulas will be useful: 



a2n = (C0S2" 



(2n- 1)!! 



d{d + 2)...{d + 2n-2) 



with n = 1, 2, . . . (obviously, ao = (1) ~ 1). From Eq. (5.24) one easily obtains: 



(sin^ d cos^" i9) 



id -I) 
{d + 2n 



■ Ol2r. 



{sm^ d cos^" d) 



id' - 1) 



{d + 2n){d + 2n + 2) 



a2n 



with n = 0, 1, 2, . . . and a2n from ( 5.24 ). 

In the notation introduced above, the integrals H2p in Eq. ( [3. If: ) are written in the form 



dk p dq 
Jm ^ Jm 'i 



sin^^z?), 35=1,2. 



Then Eq. (|5.25|) gives: 



H2 — C,i 



m-^%d-l) 
&2 



Ha — Cj 



m-'%d' - 1) 
d{d+2)e^ ■ 



(5.24) 



(5.25) 



(5.26) 



(5.27) 
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Now let us turn to the integrals h2p in Eq. (5.17). Expanding their integrands in (kq) gives: 



dq ( kq 



i+i 



(sin2Pi?cos'+iz9). 



The average (sin^^ ^? cos'"'"^ -i^) differs from zero only for odd /, and the series in Eq. ( ^.2^ ) can be written as 

oo 

h2p - -Cj In{m) (sin2p^cos2"+2^) 



(5.28) 



(5.29) 



with the integral 



Using the identity 



In{m) 



/•°° dk 


r dq 


( kq \ 


L k^+' J 




\k^ + qy 



2n+2 



(5.30) 



In{m) ^ -—Vmlnim), 'Dm = md/dm, 
Is 



which follows from the last equality in Eq. (5.3C), the integral Inijn) can be represented in the form 



In{m) 



-2e /-oo 



dk f k 



e J I fci+= \k^ + l 



2n+2 



(5.31) 



(5.32) 



where the number of integrations is reduced and the pole in e is isolated explicitly. Expanding the integrand in Eq. 
(5.32) in e and neglecting the terms of order 0{e) and higher, with the desired accuracy we obtain 



/„(m) ~ m 



-2e 



dkk^''+^ 



-Infc) = 



/l (fc2 + l)2n+2V J 4(2n+l)! 

where the quantities Bn, needed only for the three- loop calculation, can be represented as finite sums: 

(2n + l)! ^ (-1)'= [(1/2)"+'= - 1] 



6„ = ln2 



E 



n\ ^ kl{n ~ k)\{n + fc + 1)^ 



(5.33) 



(5.34) 



For the integrals /i2p in Eq. ( ^.17 ) this gives: 



/i2 = m-2-C2(d-l) ^ 



i! 2-"-2 



d(rf + 2) . . . (d + 2n + 2) 



}iA = m 



ci id' - 1) E 



i!2-"-2 



^ d(d + 2) . . . (d + 2n + 4) 



"'+S„); 



(5.35a) 
(5.35b) 



Expanding the integrand in Eq. ( 5.16| ) in 2(kq) and proceeding as above for /i2p, we obtain analogous expression 
for the integral h: 



h = m-'^ d id' ~i)y: 



i!2"-i 



d{d + 2)...{d + 2n + A) 



The 0{e ^) terms in Eqs. (5.35) can be expressed in terms of the hypergeometric series (see, e.g., [p7|) 

^, , . ^ ab a{a + l)b{b+l) 
F(a,6;c;z).l + -z+^-^-l^.- + ... 



(5.35c) 



(5.36) 



as follows: 
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h2 = Tff [-£''F{1, 1; d/2 + 2; 1/4) + (rf + 2)C2(d)] 



/14 = 
h ^ 



4d(d + 2) 

jd^ - 1) 
4d(d + 2)(d + 4) 
m-2s C| _ 1) 



2d(d + 2)(d + 4) 
and for the 0(1) terms one has: 



[-e-iF(l, 1; d/2 + 3; 1/4) + {d + 4)C4(d)] , 
[-e-^F{l, 1; d/2 + 3; 1) + (d + 4)C(d)] , 



^'^"^^ 54"(d/2 + l)...(d/2 + l + n) 

r (d]^S^ 

> ^^4"(d/2 + 2)...(d/2 + 2 + n) 



= C2(d + 2), 



C(d) = E 



(d/2 + 2)...(d/2 + 2 + n)- 



(5.37a) 
(5.37b) 
(5.37c) 

(5.38a) 
(5.38b) 
(5.38c) 



It is worth noting that the series entering into Eqs. (5.37) and ( 5.38| ) are convergent; this fact is nontrivial for h since 
1 2: 1 — 1 is the converge nce ra dius for the corresponding series. 

For F ( . .) in Eqs. (l5.37c| ) one has F{1, 1; d/2 + 3; 1) = (d + 4)/(d + 2) for any d, while the expressions for F{. . .) 
in Eqs. ( ^.37a ), ( 5.37bD simphfy only for integer d; see 27 1. In particular, for d — 2 and d = 3 one has 



d=2: F(l,l;d/2 + 2;l/4) 8[-31n(4/3) + 1] ~ 1.0956, 
F(l,l;d/2 + 3;l/4) =:6[181n(4/3)-5] ~ 1.0696, 
C2~ 0.3671, C4~ 0.2411, C ~ 0.2917; 

d = 3 : F{1,1; d/2 + 2; 1/4) 10(7r\/3 - 16/3) ~ 1.0806, 

F(l,l;d/2 + 3;l/4) = 14(-157r\/3 + 82)/5 ~ 1.0613, 
C2 =i 0.2912, C4 ~ 0.2056, C ~ 0.2412; 

and for the other integer d analogous expressions can be obtained from the recurrent relation 

3i^(l, 1; d/2 + 2; 1/4) + {d + 2)F{1, 1; d/2 + 3; l/4)/(d + 4) = 4, 



(5.39a) 



(5.39b) 



(5.40) 



valid f or all d. In Eq. (5.3£) we also included the numerical values for the coefficients C2,4 and C obtained from the 
series ( 5.3g ). 

Substituting all these expressions for the integrals into Eqs. ( 5.18 ) and ( ^.19 ) gives the final answers for the 
coefficients Ai of the diagrams Nos 22,31: 



No 22 



No 31 



A2 
A2 



1 1 1^,,, 



{p/mf^ (d^ - 1) 
16d(d + 2) 

16d^2 ' 
[ti/m)^- \ 1_ F(l,l;d/2 + 3;l/4) _ 1 



(5.41) 



4(d + 4)£ 

(d-2)(d + 3) (d- l)i^(l,l;d/2 + 2;l/4) 



24d(d+2) [e2 
[ji/mf^ (d-1) 
24d(d + 2) 



3e2 



4e 



(5.42) 



with F{. . .) from Eq. ( |5.36D and C,C2,4 from ( ^.381 ). From these expressions, using the standard general scheme, one 
obtains the contributions of the diagrams Nos 22,31 into Eq. ( 5.22| ) and then the two-loop contribution [Zp^]^ into 
the renormalization constant Zp^, presented earlier in Ref. [|oj. 
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VI. THREE-LOOP APPROXIMATION 



A. Scalarization of the three-loop diagrams 



All the needed three- loop diagrams with the symmetry coefficients are given in Fig. 1. Below we shall describe them 
in more detail and give in a compact form the complete set of relations which allow one to express the coefficients Ai 
in Eq. (4.20) for any given diagram in terms of the scalar integrals (4.41c). We recall that all the external momenta 



in the diagrams are set equal to zero. 

We begin with the normal (not factorizable) diagrams. For these, the integration momenta k, q, 1 are always 
assigned to the horizontal lines (vv) in the following order: k flows via the uppermost line, q flows via the middle line, 
and 1 flows via the lowest line. Then in order to determine the momenta for all lines in the diagram, it is sufficient 
to show (by an arrow) the chosen direction of the momentum in each line. The needed derivatives d on the lines are 
always implied (for the one- loop and two- loop diagrams, they were shown explicitly by dots in Fig. 2). T he nu merical 
indices 1,2, .. . of the upper dots are always chosen to increase from the left to the right (like in Eq. (4.17) and in 
the diagrams in Fig. 2); the positions of the letter indices a, 6, . . . in the diagram are always shown explicitly. This 
information allows one to completely restore the configuration of the momenta and the form of the index structures 
(4.39) for any diagram. 

Now we turn to the description of specific diagrams in this notation. For any diag ram, we give the directions of 
t he m ome nta, positions of the indices a,b, . . . and forms of all cofactors in Eq. (4.38) in the same form as in Eqs. 

(P-dU). 

The two-ray diagram No 23: 




No 23 

The three-ray normal diagrams: 




C = 1/2 • 1 • 1, 0, = (lPql)((q + l)Pk(q + 1)), 

= 8eieq+iek+q+i, 
/ = nl with X = k + q + 1, P = Pi. 





(6.1) 




(6.2) 



No 32 



No 33 



No 34 



No 35 



For these, one has: 



No 32 : 


C = 

4>E^ 


1 • 1 • 1, 0, = (IPql), 
4eieq+l(ek + Cq+l + Sk+q+l), 










I = 


n2 with X = k, y = q + 1, 


a -- 


= Pk(q + l), P 


= Pi [(xa) = 


No 33 : 


c = 


l-l-l, = (lPk(q + l)), 












2ei(eq + ei + eq+i)(eq + Ck+i ^ 


- Sk 


fq+l): 






I = 


n2 with X = q, y = k + 1, 


a. -- 


= Pql, P = Pi 


[(xa) = 0]; 


No 34 : 


c = 


l-l-l, 0, = (qPk(q + l)), 










4>E^ 


2ei(eq + ei + eq+i)(ei + Ck+q ^ 


- fk 


fq+l): 






I ^ 


n2 with X = 1, y = k + q. 


a. -- 


= Pql, P = Pi 


[Px = 0]; 


No 35 : 


c = 


1-1-1, <j)s = (qPkl), 










(l>E^ 


2ei(eq + ei + eq-i)(ek+q + ek+ 


1 + 


fq-l): 






I = 


n2 with x==k + l, y = q- 


1, 


a. = Pql, P = 


Pi- 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



The structure J = n2 from Eq. (4.39) is the same for all these diagrams; the corresponding quantities Ai are 
determined by the general formula (4.43) in which the specific values of the vectors x, y, z,q: and projectors P for 
any given diagram should be substituted from Eqs. (|6.3|)-(|6!6|). 
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The four-ray normal diagrams: 



No 42 

For these, one has: 

No 42 





(6.7) 



No 44 



No 45 



C=l-1-1, 03 = 1, 
(f)E^ 2ei(eq + ei + eq+i)(ek 



Cq + £1 + Ek+q+l), 



No 43 



No 44 : 



No 45 : 



/=n3 with x = k, y^q, z^l, a^P)^{q + l), P^Pql, P^Pi 

[(xa) = 0,(y/3) = 0,Pz-0]; 
C= 1/2- 1-2, 0, = (qPul), 
(/>E= 4ei(eq + ei)(eq + ei + ek-q + Ck+i), 
/ = n4 with X = q, y = 1, z = k - q, P ^ Pi, P' = Pq 

[Py = 0,P'x = 0]; 
C=l/2.(-l).2, 0, = 1, 
(/)_B= 2ei(eq + ei + ei-q)(ek + (q + ek+i + ei-q), 

/ = n3 with x = k, y = q, z = l-q, Pkl /3 = Pql, P = Pi 

[(xa) = 0,(y/3)-0,P(y + z) = 0]; 
C=l-(-l)-l, (j)s = l, 

(f)E^ 2ei(eq + ei + eq+i)(ek + ei + eq_k + Cq+i), 

/ = n3 with x = k, y = l, z = q — k, a = Pkq, /3 = Pql, P = P\ 
[(xa) = 0, ((x + z)/3) = 0, Py = 0]. 



(6.8) 



(6.9) 



(6.10) 



(6.11) 



Three diagrams have the index structure / = n3 from Eq. ( 4.39| ) and one diagram has the structure I = n4; (xa) = 
for all the diagrams with / = n3. 

Below we give the expr essio ns for the c oefficients Ai for the structures / = n3, n4, analogous to Eqs. dO^ ), ( p3| ). 

For / = n3 from Eqs. ( 1.39 ) and ( 4.40 ) using the relation (xa) = one obtains: 



12Ai = (xPx)[(ya)(z/3) + (y/3)(za)] + (yPy)(za)(x/3) + (zPz)(x/3)(ya) + 

+ 2(xPy) [(ya)(z/3) + (za)(x + y + z,(3)]+ 2(xPz) [(y/3)(za) + (ya)(x + y + z,(3)] + 
+ 2(yPz)(x/3)(z + y,a); 

72^2 = [/i{(d- l)[(ya)(z/3) + (y/3)(za)] + 2(yPa)(z/3) + 2(zPa)(y/3) + 
+ 2(yP/3)(za) + 2(zP/3)(ya) + 2(yPz)(a/3)} + 

+ C/2{(d- l)(za)(x;a) + 2(xPa)(z/3) + 2(zPa)(x/3) + 2(xP/3)(za) + 2(xPz)(a/3)| + 

+ [/3{(d - l)(ya)(x/3) + 2(xPa)(y/3) + 2(yPa)(x/3) + 2(xP/3)(ya) + 2(xPy)(a/3)} + 
+ 2(d- l)[(za)(z/3)(xy) + (ya)(y/3)(xz)] + 

+ 2[(zPa)(z/3)(xy) + (yPa)(y/3)(xz) + (xPa)(x/3)(yz) + (zP/3)(za)(xy) + (yP/3)(ya)(xz)] + 
+ 2(a/3) [(xPx)(yz) + (yPy)(xz) + (zPz)(xy)] ; 

9^3 = [/4|(d- l)(a^) +2(aP/3)|. (6.12) 
Here and below in Eq. ( |6.14[ ) 



C/i = a;2+2(xy)+2(xz), f/a = + 2(xy) + 2(yz), f/a = + 2(xz) + 2(yz), 
C/4 = a;2(yz) + y^xz) + ^^(xy) + 2(xy)(xz) + 2(xy)(yz) + 2(xz)(yz). 



(6.13) 
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For / = n4 from Eq. (4.39) using the relations Py = 0, P'x = in Eq. ( |6^ ) one obtains: 



6ii - (xPx)(yP'z) + 2(xPz)(yP'z) + (xPz)(yP'y), 
36A2 = (d- l){c/i(yP'z) + [/2(xPz) + (xPx)(yz) + (yP'y)(xz) + (z(P + P')z)(xy)} 

+ 2{j7i(yP'Pz) + C/2(xPP'z) + J73(xPP'y) + HPP' + P'P)z)(xy)}, 

9^ = C/4{(d- l)^ + 2tr(PP')} 



(6.14) 



with Ui-Ui from Eq. (|6.13| ); tr [PP') is the trace of a matrix product. 

Substituti ng th e speci fic values for the vectors and projectors for any given diagram from Eqs. ( |6.8| )-(3.11) into 
expressions (|6!T^) -(|6l4| ) gives explicit expressions for the corresponding coefficients Ai. Then, using the additional 
information from (3.j)- (|6.11 ) one obtains the desired expressions for the coefficients Ai for any diagram in the form 
of scalar integrals ( [4.41q ). Ml these simple but cumbersome technical operations are easily performed by means of a 
computer. 

Now let us turn to the factorizable diagrams Nos 46,51,61. Diagram No 46 factorizes to a product of two blocks 
Nos 21,22; diagram No 51 factorizes to a product of the blocks Nos 21,31; diagram No 61 factorizes to a product of 
three blocks No 21. Since the quantities Bi for diagrams Nos 21,22 already include their own symmetry coefficient 1/2, 
we conclude that in order to obtain the n eeded symmetry coefficient in Fig. 1 for diagrams Nos 46,51 an additional 
symmetry coefficient (like C in Eq. (5.20)) is not needed, while for diagram No 61 it should be taken to be 1/6. 

Therefore, taking into account Eqs. ( 4.18 ) and ( 4.27 ), we obtain for diagrams Nos 46,51,61: 



No 46 
No 51 
No 61 



r = 1^1234 [B1W1W2 + B2'w'^5i2\ [B[wy,Wi + B2W'^634\ , 

r = 142345 [B1W1W2 + B2w'^5i2\ [B[w:>,WiWz + B'^w'^S^iW^] , 

r = [1/6] Vi2-H5q[BiWiW2 + B2w'^5i2] [BiW^W4, + B2U?5^i\ [BiW^Wq + B2'U?5^q\ , 



(6.15a) 
(6.15b) 
(6.15c) 



where Pi, 2 are the know n coeffi cients for diagram No 21, B[ 2 in Eq. ( 6.15a ) are the analogous coefficients for diagram 
No 22, and B[ 2 in Eq. (|6.15b ) are the c oeffic ients for diagram No 31. 
Using the relations ( 4.36| ), from Eqs. ( 3.15 ) one obtains: 



r(No 46) 
r(No 51) 
r(No 61) 



kiBiB[ 
kiBiB[ 



k2{BiB'2 
k2{BiB'2 



B2B[) 
B2B[) 



k'iB2B'2, 
k3B2B' 



[1/6] [kiBf + '3k2BfB2 + 'SksBiBi + k^B^] 



(6.16) 



with the coefficients kt known from Eq. (4.36c ) fo r diagram No 46, from Eq. (4.36d) for diagram No 51, and from 
Eq. ( 4.36£| ) for diagram No 61. The relations ( |6.16 ) give the desired answers for the three-loop factorizable diagrams. 



B. Calculation of the three-loop integrals and anomalous dimensions: General scheme 



Substituting the specific values for the vectors and p rojec to rs, gi ven in ( |6.1| ), (3.2)-(6.6), ( |6.8| )-(3.11) for all normal 
three- loop diagrams, into the general formulas ( 4.42 ), ( 4.43| ), ( 3.12 )-( 6.14 ) gives explicit expressions for the quantities 
Ai in Eq. ( 4.41c| ) for any diagram in the form of polynomials in the scalar products (kq), (kl), (ql) and moduli k, q, I 
of the vectors k, q, 1. 

It is worth noting that for all the normal three-loop diagrams, the products 4>sAi have the form of sums of monomials 
in k, q, 1 of order six. In the variables "moduli-angles," each of these monomials contains the factor which is canceled 
out by the analogous factor in the energy denominator ei = i/P, present in each of these diagrams. It is also worth 
noting that these products involve the modulus k in the power 2 or less (otherwise the integrals over k would be 
divergent). 

In the three- loop diagrams, we only need the poles in e, that is, the contributions of order e~^, and e^^. The 
following general scheme is used in their calculation: 

1. The integrands are expanded in the set of scalar products (kq), (kl), (ql); the results are represented as multiple 
series in these quantities. 

2. Then, the angular averaging (. . .) is performed with respect to the directions of the vectors k, q, 1, that is, the 
following quantities are calculated: 



r„,„,„3 = ((kq)"^ (kl)"^ (ql)"^) 



(6.17) 
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with arbitrary integer exponents rii > 0. 

3. The next step is the integration over the moduh k, q, I. All the needed integrals reduce to the form: 



dk 



dq 



f^'2ni ^2n2^2n3 



dl 



(6.18) 



where the integer exponents rti > satisfy the relation ni + rt2 + 't.s =^14 + 71.5, so that the integrals contain only 
logarithmic UV divergencies for e — > 0. Only the pole parts are extracted from the integrals (6.18). 

4. The last step is the summation of the resulting series. They have the forms of double infinite series with the 
coefficients given by n-fold finite sums [n < 5); the number of terms in the latter increases rapidly with the order 
of the coefficient. This summation is the only operation which cannot be performed exactly (analytically) and is 
therefore the only source of errors in numerical coefficients in expressions like Eqs. (1.12). 

Of course, the straightforward but cumbersome operations listed above have been performed with the aid of a 
computer. 

The first step, the expansion in the scalar products (kq), (kl), (ql), contains some conceptual subtleties, and we 
shall discuss it in more detail. The problem is that the plain expansion of the integrands in the powers of the scalar 
products in some cases leads to divergent series. 

The non-polynomial dependence of the int egr an ds o n the co sine s of the angles appears only from the energy factors 
(/)^^, whose explicit forms are given in Eqs. ( |6.lD , (^3|)-(6.f:), (3.8)-(6.11). For all the diagrams, the nontrivial factors 
in <j)E have the forms + 1^ + const (ql) and k^ + q^ +lr+ some linear combination of all scalar products, and in 
the expansions in powers of the scalar products, the denominators will contain powers of the quantities (g^ + P) and 
{k^ + q'^ + P). Such expansions converge for all cofactors in 0^^, which do not contain the "energy" Ck+q+i with the 
sum of all three momenta. The cofactors with ek+q+i require special consideration. They are present in diagrams Nos 
23,32,33,34,42 and are proportional to the following factors: 



No 23 
No 32 
No 33 
No 34 
No 42 



6^^ CX [ek + Eq+l + Ck+q+l]"^ cc [Q + S - 
OC [eq + Ek+l + Ck+q+l]"^ CX [Q + S" - 



1)^ (X [ei + ek+q 



Ck+q+l] ^ OC [0 - 



s 



(kq)]-\ 



6^^ CX [ek + Eq + ei + Ck+q+l] ^cc[Q + S] \ 



(6.19) 



where Q = k^ + q^ + P and S = (kq) + (kl) + (ql). 

From the obvious identities (k + q + 1)^ > 0, (k — q)^ 
it follows that 



(k-l)2 + (q-l)2 > 0, |(kq)| <kq, |(kl)| < fc^ and |(ql)| < ql 



Q/2<S<Q, |(kq)|<Q/2, |(kl)|<Q/2, |(ql)|<Q/2. 



(6.20) 



It then follows that the factor ( 6.19| ) in diagram No 42 can be expanded in the powers of the ratio S/Q with 
\Si/Q\ < 1, while for No 23 this is impossible: in the integration region, there is a subregion of the same dimension 
(namely, 3 x c?) in which 2151 > Q. This difficulty can be circumvented by means of the shift of the point around 
which the expansion is made in diagram No 23: 



[Q + 2S]-^ ^ [3Q/2 + (45 - Q)/2]-i = [2/3Q] ^[(1 - 45/Q)/3]' 



(6.21) 



1 



The c onvergence of the series in (6.21) is ensured by the inequality 
( |6.20 ). The equality takes place only in the 2 x c?-dimensional subregion k 



- 4S'/Q)/3| < 1, which follows from Eq. 
■ q-|- 1 = 0, which has zero measure in the 
3 X d-dimensional integration region and does not spoil the conve rgence of the corresponding series for the integral 
(similar consider ations ensure the convergence of the series (5.35c) for the two-loop integral h). 
For the factor (6.19) in diagram No 32 the following shift can be used: 



[Q + S+ (ql)]-i = {5Q/4 +[S+ (ql) - Q/4]}-' 

with the subsequent expansion in the powers of the ratio [S + (ql) — (5/4]/ (5(3/4), whose modulus is < 1 owing to 
inequalities (6.20). Similar expansions (with obvious modifications) can be written for diagrams Nos 33,34. 

Such infinite series contain additional finite sums originated from the powers of the expansion parameters, which 
have the forms "constant + linear combination of the scalar products." It is clear that this summation should be 
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performed in the first place, before all the other summations: only this order of summations ensures the convergence 

of the series. 

Calculation of the angular (6.17) and momentum ( 3.1^ ) integrals is a separate task; it is discussed in the next two 
Subsections. The results for the coefficients Ai for all the normal three-loop diagrams are presented in the Appendix. 
The e^^ contributions in all diagrams have been found analytically for general space dimensionality d (in order to 
calculate them, it is sufficient to neglect all the scalar products in the energy denominators <j>E), while the and 
parts have been found for the physical cases d — 2,3 and in the large d limit; the results have been presented in 

Ref. (of course, this calculation can be made for any fixed given value of d). 

Now, using the standard scheme (see Sec. IV), one finds the three-loop contribution [^i^^jg in the expansion (4.7) 
for the re norm alization constant Z^ ^. T he ano mal ous dimension 7^? = jni is given by the relation — —(idu \nZp^; 
see Eq. (3.4). Substituting Eqs. (2.17) and (4.7) into the last relation expresses the anomalous dimension in the 



coefficients [Zp^]^^. Within our accuracy one obtains 

Mu)=[e^u{d-l)/2d]{[Zp']^+2[Zp']^-[Zp']l+3[Zp\-3[Zp']^[Zp\+[Zp']l}. 



(6.22) 



Substituting the known expressions for [Zp^]^ into Eq. (|6!22|) gives the anomalous dimension 7f(u) to order u 
Since the quantity 'jp is UV finite, that is, finite at e = 0, the pole parts must cancel each other in Eq. (3.22) 



This implies some exact relations between the senior poles in the quantities [Zp^]/^ with k > 2 (e^^ in [Zp^]^ and 
in [^^j^^jg) and the parts of the previous orders in u. Such relations provide an additional possibility 
to control the absence of calculational errors. In fact, the knowle dge o f senior poles (e"'^ with fc > 2) is needed only 



£-2, £-3 



to check this cancellation; the nonvanishing contributions in Eq. (3.25) are completely determined by the e terms 



in the quantities [-^^i^^Jj,- In the MS scheme, the anomalous dimension 7_f(u) appears independent of e; this is a 
consequence of the rela tion Zp — 1+ only poles in e . 

Finally, the relat ion (|3.4| ) with from Eq. ( 2.18| ) gives the 0{e^) results for the critical dimensions presented in 
Ref. [H and Eqs. (pTlS). 



C. Angular integrations in the three-loop diagrams 



The two- loop integrals (4.41b) involve only two vectors k and q and o ne an gle ?? between them, so that the pr ocedure 
of the angular averaging reduces there to the only standard formula ( 5.24 ). The three- loop inte grals ( 4.41c ) involve 
three vectors k, q and 1 and three angles between them, so that the calculation of the quantities ( 6.17 ) is not all that 
simple. Below we present t he re sults of this calculation. 

Obviously, the quantity ( 6.17 ) differs from zero only if all the three numbers 71.1.2,3 are simultaneously even or odd. 
It is also clear that the quantity (3.17) is symmetrical with respect to any permutation of the exponents ?ii,2,3, and 
with no loss of generality it can be assumed that n2 is the minimal exponent: 



712 = min{7ii, 7i2, 713} 



(6.23) 



(the notation in the formulas below is consistent with the assumption ( 6.23 )). 

The straightforward calculation (first, the averaging over the direction of one momentum, say k, and then the 
avera ging over the angle between the two remaining vectors q and 1) leads to the following result for the quantity 
(pTh: 



T — l,ni+n2 „ni+n3 /n2+n3 rp 

where the moduli of the vectors are isolated explicitly, and 

^711712^3 — Q^ni +n2 '-^n2-|-n3 ^^(^1 ; 7l2 ) ^nin2n3 

with coefficients a2n from Eq. ( ^.24 ) and 

2"2 nil [{ni + n2)/2]\ _ (tii - n2 + 2){ni - 712 + 4) . . . (ni 712) 



K{ni,n2) = 



(ni+n2)![(7ii-n2)/2]! 



{ni + l){ni+2)...{ni+n2) 



(6.24) 



(6.25) 



(6.26) 



T, 



712/2 



nin2n3 ~ / , 
k=0 



Yk 



I ^''712! [(tii - 7i2)/2]!a„2+„3_2fe 



k\ (712 - 2fc)! [(ni - n2)/2 + k]\a. 



(6.27) 



"2 +"3 
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We rec all th at the numbers n.i_2,3 in these e xpres sions have the same parity, and that n,2 is the minimal one according 
to Eq. ( 6.23 ). The upper limit in the sum ( 6.27| ) for odd 712 is understood as the in teger part of 712/2. 

The following special values and recurrent relations for the quantity K in Eq. ( 6.26| ) are useful in the computer 
calculations: 

A'(m,0) = K{m, 1) = 1, K{m,m) = ml/{2m- 1)!!, K{m+ l,m + 1)/ K{m,m) = (m + l)/(2m + 1), 



K{m + 2,n)/K{m,n) = {m + l)(m + 2) / (m - n + 2)(rn + n + 1). 
For the terms Yk in the sum ( |6.27 ) one has: 

(na - 2fc)(n2 - 2fc - + ria + - 2fc - 2) 



Yo = 1, Yk+i/Yk 



2{k + l)(ni -n2 + 2k + 2)(n2 + 773 - 2A: - 1) 
Using these relations, the quantities (6.26), ( |6.27| ) are easily calculated by a computer. 



(6.28) 



(6.29) 



D. Modular integrations in the three-loop diagrams 



Let us turn to the calculation of the th ree-lo op integrals (|6.18|) over the moduli fc, q, I, whic h aris e as coefficients in 
the expansion of the quantities Ai in Eq. (4.41c) in the scalar products (kq), (kl), (ql); se e Sec . VI B, For the complete 
three-loop calculation, it is sufficient to know some special cases of the general integral ( 6.18 ), which we denote /1-/9 
in what follows. Below we give the explicit answers for these intergals, with the precise specification of the indices 
711-715, and then turn to the derivation. Only the pole parts £~^, and of these integrals are given, which is 
sufhcient for the calculation of the quantities Ai within our accuracy. We shall use the notation 



for any function F = F{k, q, I). The integrals /1-/9 are as follows: 

r P"i^2n2;2n3 ] (77i - 1) ! (?^2 - 1) ! (7^3 - 1) ! ("4 " - 1) ! 

|(fc2+g2_^P)«4(g2+^2)„5 J - 12£ (774 " 1) ! (7l2 + «3 " 1) ! 

with 77i + 772 + 773 = 714 -I- 715 , 774 > ni , «i,2,3,4 > and 715 > 0; 

g2"2;2«3 I - l)!(7i3 - 1)! 



(6.30) 



(6.31) 



(fc2+g2 + p)n4(g2_^;2)n5 



12(772 + 773 - 1)! 



^ 2 /"2-(-n3-l 2 "4 — 1 2 ^ n2-l ^ ^ 113 — 1 ^ 



^ k ^ k 

k=l k=l 



2 ^ /c 

fe=i 



2 ^ fc 

fe=i / 



(6.32) 



with 772 + 773 = 774 + 715 , 772,3,4 > and 775 > (in Eq. ( 6.32 ) and all the formulas below, any sum with the upper 
limit lesser than the lower one is understood as equal to zero); 



fe2"ig2n2 



(fc2 + q^ + 
with 77i + 772 = 774 and 711,2 > 0; 

^2nig2n2 



(fc2-Fg2+;2)n4(q2+;2) 

with 77i + 772 = 714 + 1, 77i > and 772 > 1; 

fc2 f ql 



(771 - 1)1(772 - 1)! 

12(774 - 1)! 



(771 - 1)!(772 - 2)! 

12£ (774- 1)! 



p2 Op 7, ^ 2^ 7. 



£2 2e \ ^ k 

\ k=l 



1 1 



k=l 



(6.33) 



2e I (772 - 1) 



ni-l ^ n2-2 ^ 

^k^^k 

fe=l k=l / 



(6.34) 



(fc2 + g2+p) \qi + l 



2n^ 



12(277 - 1)! 



2n-l 



k—n 



(6.35) 
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with n > and Bn from Eq. (|5.34|) ; 



fc2g2(n+l);2n 



(fc2 +q2 +;2)(q2 _^;2)2n+l 



(6.36) 



/7 = / 



A:^ + g2 + p 



1 

3?' 



(6.37) 



2„2 



(fc2+92+p)(52_^p) 



6e-^ ' 



(6.38) 



(fc2+g2+;2)(^2+p) 



1 

3? 



1 

T2s' 



(6.39) 



Now let us turn to the derivation of expressions ( 3.31 )-( 6.39 ). The integrals /i,2 can be obtained from the generation 
function 



i?(a, 6; m) 



dk 



dp 



-3e r 



We shall also need the first terms of the e expansion of the integral 



\ + - Hb/a) + 0(1) 



(6.40) 



JiP,q) = / dip {cos ip) 
Jo 



2p-l-£ 



(sin ip) 



l + e 



2T{p + q) 

(p-l)!(g-l)! 
2{p + q-iy. 



i'iv + q) - ^^{p) - ^VX?) 



r(p-£/2)r((7-£/2) 

2rip + q-e) 

9-1 



E 



k=l k=l k=l / 



(6.41) 



where V{z) is the Euler Gamma function and "^{z) = c?lnr(z)/dz. 



The calculation of the quantity i? = i?(a, 6; m) 



i?(a, b; 1) is s imilar to the derivation of Eq. (5.33). The 



operation 'Dm = md/dm is applied to the double integral in Eq. (6.4C), which reduces it to the sum of two single 
integrals and explicitly isolates the pole factor £~^: 



-3e 



R 



3e 



-T>,„R 



3e 



dk- 



dp 



k^+^{ak'^ + b) p^+^'' {a + bp^) 



Now one can set £ = in the integral over k (it remains finite at £ = 0). The pole in e in the integral over p comes 
from large p and it can be isolated explicitly: 



dp 



1 



^1+2-^ (a + 6p2) 5 



dp 

nl + 2£ 



dp 



1 



1 



b J I (a + 6p2) 2be 



dp 



p{a + bp^) 



+ 0(£) 



which immediately leads to the answer ( 6.40| ). 

In the integral Ii the pole contribution comes from the part of the integration region where all three momenta 
k, q, I simultaneously tend to infinity, and in I2 there is also a pole contribution coming from the subregion where q, I 
simultaneously tend to infinity at fixed finite k. This means that the pole parts of the integrals /i_2 are not affected 
if one changes to the polar coordinates p, ip {q — pcosip, I — p sin 1^9) with the integration region < < 7r/2, 
1 < p < 00, that is, the following replacement is performed: 

POO /'7r/2 POO 

dq dl . . . ^ dp pdp . . . . 

1 Ji Jo Ji 

The integration regions in these two expressions are not identical, but the values of integrals differ only by an 
unessential contribution, finite at £ ^ 0. For Ii this gives: 
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h 



tt/2 



dip (cos Lp) 



2n2 



-1-e 



(sin if) 



2n, 



dk 



dp fc2"ip2«4-2ni 



nl + 2£ 



(fc2 + p2) 



2VJ4 



This integral can be expressed using J(p, q) from Eq. (6.41) and the derivative of the generating function R{a, b; m) 
from Eq. (3.40) as foUows: 



h = [(^4-1)!]-^ J(n2,n3)(-aa)"U-<9,)^ 



-^R{a,b;l)\a=b=i 



(6.42) 



with da — d/da, db — d/db, whi ch immediately leads to the answer ( 6.31 ). Note that the e ^ contribution in R 



depends only on b (see Eq. (6.40)) and therefore vanishes in the above expression owing to the d iffere ntiation with 
respect to a (we recall that ni > for Ii). Although the integral I2 is formally given by expression ( 6.42 ) with ni ~ 0, 



h = [(n4-l)rij(n2,n3)(-5b)"^-ii?(l, 6; 1)1,^1 



(6.43) 



the part of the function (|6.40|) survives in Eq. ( 6.43 ) owing to the absense of the derivative da, so that the 0{e) 
contribution s sho ul d be t aken into account in the integral (3.41). This explains the sharp difference between the final 

expressions ( |6.31 ), ( 6.32) for the integrals Ii and l2- 

The relations (3.33)-( 6.36| ) follow easily from (3.31), ( |6.32[ ). The result (| 6.33| ) for I3 is obtained fror a En . ( 6.32 ) 
at ^5 = using the symmetries of the integrand. In its turn, the result ( |6.34| ) for J4 follows from Eqs. ( 3.31 ), ( 6.32 ) 
i f the obvious identity (?2"V(9^ + l^) = g^^"""^) - q2("2-i)p/(-^2 ^ substituted on its left-hand side. The result 



(3.35) for J5 is obtained using the substitution fc^/(fc^ + q'^ + l'^) = 1 — (q^ + l'^)/(k'^ + q^ + l'^)- i n the fir st ter m, the 
integral over k is trivial and the remaining integrations over g, I are performed using Eqs. (Isisq) , ( 15^3^ ), (I^H), while 
the second term reduces to l2- The integral /g in Eq. ( |6.3(^ ) reduces to /s after the symmctrization of the integrand 

with respect t o g, I. 

The results (6.37), ( |6.38 ) follow from the obvious relation /{I} 



e in Eq. (6.30) after the symmctrization of the 



integrands with respect to fc, 5, 1. The last relation (|6.39D is obtained using the substitution q^ / {q'^+P) — l—l^/[q^+l^) 
the first term coincides with /y and the second one is given by I2 with n2 — ~ = nc, = 1. 



E. Anomalous exponents at large d 



In th is Se c tion , we shall briefiy discuss the behavior of the coefficients A^^'' in the expansion (1.10) for d 00. The 
model (1.1 )-( fL3| ) has no finite upper critical dimension, above which the anomalous scaling vanishes. It disappears 
for infinite d pSf, but reveals itself already in the 0{l/d) approximation B . This fact confirms the relevance of the 
large d expansions for the issue of anomalous scaling; see also Refs. [^,^ for the discussion of the Navier-Stokes 

problem. 

It follows from Eqs. (5.24) that for d 00, the angular averages a2n = (cos^""!?) behave as a2n d^" (each 
additional factor cos^ introduces additional smallness in 1/d). Then from Eqs. (6.24)-( p.27] ) for the averages (6.17) 
one obtains T„i„,„3 cx d-("i+"2+"3)/2 foj. „^ ^,3 even and T„i„2„3 cx d-("i+"2+n3+i77^or ^1,2,3 odd. This means that, 
in order to find the behavior of the coefficients Ai for large d to any given fini te orde r in 1/d, one needs to take into 
account only /tmte numb er of terms in the expansion of the integrands in Eq. (4.37b) in the scalar product (kq) and 
integrands in Eq. ( 4.37c ) in the set of scalar products (kq), (kl), (ql). 

The d dependence of the coefficient A^"';-' in Eq. ( 1.10| ) is known from Eq. ( 1.11 ), while the quantities A^^^;'' 
be found as series in l/d to any given order as explained above. For general n and I to order l/d^ we have obtained 



I : ^ni can 



A„, = e[ - n{n - 2)(1 - 2/d)/2d + {1/2){1 - 2/d + l/d + 2/d^)] + 3e^{n - 2){n ~ l)/4:d^ + 
+ £^(71- 0[l-74988(n- 2) - 0.624916 Zj/d^ + 0{e^). 



(6.44) 



Note that the and contributions d ecay for d ^ 00 faster than l/d in agreement with the 0{l/d) result obtained 
in Ref. for A„o- Moreover, from Eq. ( 6.44 ) it follows that the leading 0{l/d^) terms in these contributions vanish 
for n — I, so that the decay at d — s- cx) becomes even faster. For n — I we have obtained 



A, 



en 



/2 + n{n - 1) {e /{d -l){d + 2)- [l + {2n - 7)/d] /d^ - e^{3n ~ 8)/2d^} + 0{e^) 



(6.45) 



with the accuracy of 0(1/ d'^). 
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VII. CONVERGENCE OF THE e EXPANSION, INVERSE e EXPANSION, AND COMPARISON WITH 

NONPERTURBATIVE RESULTS 



An important issue which can be discussed on the example of the rapid-change model is that of the nature and 
convergence of e expansions in models of turbulence and the possibility of their extrapolation to finite values of 
e ~ 1. The knowledge of the three terms of the e expansion in model (|n|)-(^) allows one to discuss its convergence 
properties and to obtain improved predictions for finite e in reasonable agreement with the existing nonperturbative 
results: analytical solution of the zero- mode equations for n = 2 Q, numerical solutions for n = 3 and numerical 
experiments for n = 4 ||2^,|2^ and n — Q psf . 

In Figs. 3a and 3&, we show the anomalous dimension 7 = 722 = 722('"*) for d = 3 (a) and 2 (&) in the 0(e), 0{e^) 
and 0(£^) approximations (from above to below); the latter two obtained as simple sums of two and three terms of 
the e expansion, respectively. The dashed line corresponds to the exact solution by Ref. j|] ; see also Ref. |^ for the 
special cases d = 2 and 3. Analogous diagrams for the cases n = 3, Z = 1 and n = 4, Z = can be found in Ref. p5| |. 

All these figures show that the agreement between the e expansion and nonperturbative results for small e improves 
when the higher-order terms are taken into account, but the deviation becomes remarkable for e ~ 1 and decreasing d. 
Furthermore, the convergence of the e series appears more irregular for d = 2, while the forms of the nonperturbative 
results are not much affected by the choice of d. 

Such behavior can be naturally understood on the example of the exact analytical result for 7 = 722 @ , which can 
be written in the form 

27= -(d + 2 + e) + V(e + £+)(£ + £-), e± = (^d^ + d+2± ^/8d{d + 1) ) /{d - 1) (7.1) 
with e+e_ = {d + 2)^. It is useful to rewrite the smaller quantity e_ in the form 

e_ = (d-l)|l+(l + 3d+v/Mt^TT)) (7.2) 

From these expressions it follows that the corresponding e expansion has the finite radius of convergence e_, ranging 
from to 00 when d varies from 1 to 00 (in particular, e_ ~ 1.1, e+ ~ 14.5 for d = 2 and e_ ~ 2.1, e+ ~ 11.9 for 
d = 3). Hence, the naive summation of the e expansio n fo r 7 works only in the interval e < e_, which decreases 
almost linearly with (d—l). Since the singularity in Eq. ( [7.l| ) occurs for negative e, it affects strongly the convergence 
of the e expansion but is not "visible" in the form of the exact curve for positive e, in contrast with the singularities 
occurring at e = 2 in higher-order critical dimensions [[7||^ j2^ , |2^ . 

Therefore, in order to recover the behavior of 7 from its e series for larger e, it is necessary to isolate explicitly the 



singularity at e = — e_ in Eq. (7.1), thus changing to a kind of improved e expansion, whose radius of convergence is 
determined by a more distant singularity. This can be done , fo r example, by introducing the new expansion parameter 
X = + £- — ^/e^, that is, e — x"^ + 2x^/£^. Then Eq. ( [7.l| ) can be written in the form 



27 = - (d -I- 2 + 2x./£Z + x^) + {x + ^/?r) + 2x^/eZ + e+ . (7.3) 

The convergence radius of the expansion in x is found from the equation -t- 2x^f£^ -I- e+ = with the solutions 
x± = — y^eZ ± 1^/6+ — £- and is therefore equal to \x±\ — y^e^. For positive x and e, this corresponds to the 
convergence for < a; < y/e+ or, equivalently, < e < e+ + 2y/e-e+ = £+ + 2{d + 2). 

The improvement of the convergence is illustrated by Fig. 3 c, where the exact exponent 7 for d = 2 is shown as a 
function of e along with its first (x), second {x and x'^) and third {x, x^ and x^) orders of the improved x expansion, in 
which the variable x is also expressed as a function of e. One can see that the convergence of the x expansion appears 
more regular and that its third-order approximation is hardly distinguishable from the exact result for all < e < 2 
(for d — 3, the agreement is even better and for this reason is not shown). It should be stressed that crucial for the 
improvement was not the existence of the exact solution or explicit form of the substitution x{e) but the knowledge 
of the character and location of the singularity which determines the convergence properties of the plain e expansion. 

The difference with the models of critical phenomena, where e series are always asymptotical, can be traced back 
to the fact that in the rapid-change models, there is no factorial growth of the number of diagrams in higher orders of 
the perturbation theory. The divergence for c? — > 1 is naturally explained by the fact that the transverse vector field 



does not exist in one dimension. We also recall that the RG fixed point diverges at c? = 1; see Eq. (2.18), so that the 
coefficients of the e series diverge for all dimensions A„; . 

It is then natural to assume that the e series for all A„/ also have finite radii of convergence with the behavior similar 



to that of e_ in Eq. (7.2). Therefore, in order to improve their convergence and to obtain reasonable predictions for 
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finite values of e, one should augment plain e expansions by the information about the character of the singularities 
and their location in the complex e plane. Such information can be extracted from the asymptotical behavior of the 

(k) A K 

coefficients A^; in Eq. (1.10) at large k. To our knowledge, the large k behavior of the e series remains an open 
problem for any dynamical model; the instanton analysis developed in Refs. ]3^ ] has mostly been concentrated on 
the behavior of the exponents in the limit n — > oo. One can hope that the implementation of the instanton calculus 
within the RG framework will give the solution of this important problem. 

It turns out, however, that certain elementary considerations allow one to improve the convergence of the e series 
and, at the same time, to achieve a better agreement with the nonperturbative results. Let us explain the idea on 
the example of the exact solution (7.1). We express e as a function of the exponent 7 using Eq. (7.1) and expand the 
right-hand side of the resulting exact relation 



e = 7 (d - l){d + 2 + 7)/[2 - 7(d - 1) 



(7.4) 



in 7; this gives the "inverse 7 expansion." It is easy to see that, for e ~ 1 and physical dimensions d — 2 and 3, the 
respective values of 7 lie within the region of convergence of the inverted series: 7 < 2/{d— 1). The improvement of 
the convergence is also seen from Fig. 3(i, where the exponent 7 = 722 for d = 2 is shown as a function of e along with 
the first (7), second (7 and 7^) and third (7, 7^ and 7^) orders of the 7 expansion, expressed in the origin al v ariable 
e: the approximate curves approach the exact curve from the same side and represent the exact result (7.1) much 
better than the corresponding approximations of the plain e expansion. The improvement is even better for d = 3. 

A simple explanation of the improvement foll ows. The convergence radius for the direct e series is determined by 
the singularity in the right-hand si de o f Eq. (7.1), closest to the origin, that is, Ec — — £_. This square-root singularity 
disappears in the inverse relation (7^), that is, the dependence of e on 7 in the vicinity of the corresponding point 
27c = —d — 2 — £c — —d — 2 + E- becomes analytic. This would also happen for any singularity of the form 
7 — 7c oc (e — EcY^'' with any integer fc > 0. 

We assume that these features are also typical to the higher-order dimensions A„; and construct the corresponding 
7 expansions. It is important here that 7 = 7^; = 0{e), so that there is a one-to-one correspondence between these 
two expansions and three terms of the 7 expansion can be immediately obtained from Eqs. (1.7), ( |1.12D and (1.13). 
This simple procedure leads to remarkable improvement of the convergence and, at the same time, the agreement 
with the numerical results, as is easily seen from Figs. 4 and 5. 

In Fig. 4, we show the anomalous dimension 7 = 73]^ for d = 3 (left) and d = 2 (right): the 0{e) approximation, 
the third-order approximation of the inverse 7 expansion, and the third-order approximation of the plain e expansion 
(from above to below); the dashed lines represent the exact numerical solution by Refs. ||^ (7 = A — 3 in the notation 

ofi).^ 

In Figs. 5a and 5b, we show the quantities A„ = A„o in three dimensions for n — 6 (a) and n = 4 (b): the 0{e) 
slope, the third-order approximation of the 7 expansion, and the third-order approximation of the plain e expansion 
(from above to below) ; the dots connected by dashed lines represent the results of the numerical simulations by Refs. 
[H (n = 6) and |2|,|^ (n = 4). 

In all these cases, the improvement in the agreement with nonperturbative results is obvious. It should be empha- 
sized, however, that even the plain e expansion captures some qualitative features of the dimensions A„ established 
in the numerical simulations p^-^25|]: the quantity |A„| increases with e, achieves a maximum at some point inside 
the interval < e < 2 and then decreases to zero; the height of the maximum increases and its position moves to the 
left as n grows from 4 to 6 or d decreases from 3 to 2 {(^2n — 'nC,2 = ^2n in the notation of Refs. ||2^-|25|). 

It is no surprise, of course, that the disagreement between the perturbative and nonperturbative results becomes 
rather strong for e > 1; this can be explained, for example, by the effect of the singularity at e = 2 in the exact 
solutions [[7||q, |23| , p^ and by insufficient number of the known terms in the e and 7 series. For the case n = 4, d = 3, 
the situation can be improved by an interpolation formula which takes into account the first terms of the e expansion 
along with the asymptotical behavior of the dimension A4 in the vicinity of the opposite edge e = 2, known from the 



numerical simulation 



[0.06 {2-e) + 1.13 (2 - ef/"^] . In particular, one can choose: 



A4 = - 



ci£ + C2e^ + cse-' + . . . + + 



0.06 (2 -e) + 1.13 (2 -e) 



3/2 



(7.5) 



The first coefficients ci-Cp are determined by the requirement that the expansion in e of the right-hand side of Eq. 



(7.5) reproduce the first p terms of the e expansion for A4, known from the RG (therefore, in practice one can only 
take p < 3). The values of these coefficients, once determined, will not change if one takes a larger value of p. The 
remaining coefficients Ck with k > p should be chosen to reproduce the correct behavior at e — > 2. The simplest 
possibility is to set Ck = for all k > p + 1; then the last remaining coefficient Cp+i is unambiguously determined by 
the relation (2ci -I- 4c2 -I- 8C3 -I- . . . -|- 2P+^Cp+i) = 1 [the 0{2 — e) terms in the expression in the first square brackets 
produce only corrections of order (2 — e)^ to the behavior of the right-hand of Eq. (7.5) side at e 2 and thus 
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should be neglected]. This procedure gives c\ = 0.241, C2 = 0.129 for p = 1 (upper curve in Fig. 5c) and c\ = 0.241, 
C2 = 0.168, C3 = —0.225, C4 = 0.103 for p = 3 (lower curve). One can see that the inclusion of the higher orders 
of the £ expansion improves remarkably the agreement with the "experimental" results of Refs. [^,|^, shown in 
Fig. 5c by the dots connected by dashed lines. It is also worth noting that the value of the coefficient ci for p — 1 
(determined by the behavior at e — > 2) appears rather close to its value for p = 3 (determined by the e expansion), 
which demonstrates the robustness of the results obtained by the above procedure. 




(a) (b) (c) (d) 

FIG. 3. Anomalous dimension 7 = 7I2 for d = 3 (a) and d = 2 {b,c,d). Dashed line: exact solution by Refs. ^Jsi). Solid lines 
in Figs. 3a and 36: first, third and second approximations of the plain e expansion (from above to below). Solid lines in Fig. 3c: 
first, third and second approximations of the improved x expansion (from above to below; the third approximation is practically 
indistinguishable from the exact solution for all < e < 2). Solid lines in Fig. 3d: first, second and third approximations of the 
inverse 7 expansion (from above to below). 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



FIG. 4. Anomalous dimension 7 = 7J1 for d — 3 (left) and d = 2 (right): the 0{e) approximation, the third-order 
approximation of the inverse 7 expansion^and the third-order approximation of the plain e expansion (from above to below). 
Dashed line: numerical solution by Refs. |7|. 




FIG. 5. Critical dimensions A„ for d = 3: n = 6 (a) and n = 4 (6,c). Dots connected by dashed lines: numerical simulations 
by Refs. (n = 6) and |23[Q (n = 4). Solid lines in Figs. 5a, b: the 0(e) slope, third-order approximation of the 7 expansion, 
and third-order approximation of the plain e expansion (from above to below). Solid lines in Fig. 5c: approximations for n — 4, 
obtained using the interpolation formula (7.5) with p — 1 (upper curve) and p = 3 (lower curve). 
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VIII. CONCLUSION 



To conclude with, we have studied the inertial-range anomalous scaling of a passive scalar quantity a dvec te d by the 
Gaussian velocity field, white in time and self-similar in space. The corresponding stochastic problem (l.l)-( |l.3| ) can 
be reformulated as a field theoretic model (2.1), which allows one to identify the anomalous exponents with the critical 
dimensions of certain scalar and tensor composite operators built of the scalar gradients and to calculate them within 
the RG and OPE approach in the form of a regular perturbation expansion, similar to the well-known e expansion in 
the RG theory of critical behavior. 

Earlier, the anomalous exponents were presented to orders |]l0|- |l2| , |l^ and the main goal of the present 

paper has been the detailed explanation of the corresponding calculational techniques and derivation of the three-loop 
result, including the anisotropic sectors. Owing to the comparative universality of the RG and OPE formalism, these 
techniques can be applied to other models of dynamical critical phenomena and systems far from equilibrium: passive 
advection by the non-Gaussian velocities with finite correlation time, stochastic Navier-Stokes equation and so on; 
see Refs. 

Another scope of the paper has been the discussion of the convergence properties of the e expansion and the 
possibility of its extrapolation to finite values of £. It was shown that the knowledge of three terms allows one to 
obtain reasonable predictions for finite £ ~ 1; even the plain e expansion captures some subtle qualitative features of 
the anomalous exponents established in numerical experiments. 

We believe that the framework of the renormalization group and operator product expansion, the concept of 
dangerous composite operators and the e expansion will become the necessary elements of the appearing theory of 
the anomalous scaling in fully developed turbulence. 
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APPENDIX A: COEFFICIENTS Ai FOR THE THREE-LOOP DIAGRAMS 



Below we give the pole parts of the coefficients Ai from Eq. ( 1.2C| ) for all normal (not factorizable) three-loop 
diagrams. We use the following notation: 



A- = w 



Coefficients A^^'"^^ are given in Tables 1 and 2 for d = 2 and o? = 3, respectively; coefficients ai{d) are given in Table 3 
for general d. 







TABLE L 


Coefficients 


for d = 2. 






No 


Af^ ■ 10^ 


A« ■ 10-' 


Af . 10^ 


A'^) • 10^ 


Af ' • 10^ 


A<i) • 10^ 


23 
32 
33 
34 
35 
42 
43 
44 
45 


0.23875 
-0.56396 

1.11953 
-1.05470 
-0.41468 
-0.00053 

0.02902 
-0.03113 
-0.26802 


-0.06103 

0.32532 
-0.49834 
-0.12495 

0.30671 
-0.11026 

0.00249 
-0.01477 

0.26530 




0.49056 
0.93382 
0.2434 

-0.63626 
0.06564 
0.03161 
0.03350 

-0.02352 




-0.56645 
-1.04911 
-0.5770 

0.22330 
-0.16222 
-0.00715 
-0.00286 

0.03554 


-0.54061 
-0.33590 
-0.19626 
0.88599 


0.32914 
0.07208 
0.21410 
-0.88934 
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TABLE II. Coefficients '^^ for d = 3. 



No Af^ ■ 10^ A^^^ ■ 10^ ■ 103 A^^^ ■ 103 Af^ ■ 10^ A^^^ • 10^ 

23 1.03775 -0.14683 — — 

32 -0.52289 0.39815 0.21884 -0.50345 — — 

33 1.57028 0.00156 1.56822 -0.76955 — — 

34 -1.64375 -0.82315 -0.43185 -1.46450 — — 

35 -0.43289 0.51544 -1.07015 0.37916 — — 

42 -0.07198 -0.08409 0.02952 -0.22968 -0.83453 0.37471 

43 0.02808 0.00152 0.03587 -0.01186 -0.62532 0.07298 

44 -0.07726 0.01392 0.00504 0.02341 -0.36121 0.37991 

45 -0.33881 0.33509 -0.10389 0.11756 0.94431 -0.95923 



TABLE III. Coefficients ai{d) for general d. 

No ai{d) 0-2 (rf) asid) 

23 9(d+l)(2d^ + 2d-5) 27(d - l^id + 2)'^ = 

32 6(d + l)(d2 +rf-3) 2(d* + 3^3 - 3d^ - 7d + 10) — 

33 -12(d + l) -4(d2 + 2d-2) — 

34 3(d + 1) (3d^ + 3d - 8) 3d* + 8d^ - 13d^ - 22d + 32 — 

35 -12(d + l) -4(d^ + 2d-2) — 

42 3(d+l)(d^+d-8)/4 (d* + 6d3 + 5d^ - 16d + 32)/8 -4(d^ + 2d - 2) 

43 -6(d + 1) (d^ + d^ - 4d + 8) /2 d* + 4d^ + d^ - 6d + 8 

44 -6(d+l) (d^ + d^ -4d + 8)/2 d"* + 4d3 + d^ - 6d + 8 

45 3(d+l)(d^+d-8)/4 (d'' + 6d3 + 5d^ - 16d + 32)/8 -4(d^ + 2d - 2) 
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